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PREFACE 


At the present time, when the real possibilities of conquering cosmic 
space have been disclosed by the world-renowned achievements of Soviet 
science and engineering, there is great interest in the dynamic problems 
of near-earth and cosmic flights. It is expedient to start the study of 
complex space flights with the study of Keplerian motions, the simplest 
kinds of motions in cosmic space. Among a large number of readers 
interested in the problems of flight in outer space, the need has been felt 
for a small book on the theory of Keplerian motions, in which all kinds of 
Keplerian motions and their characteristics that are of interest in con- 
nection with space flights would be considered in a compressed but 
sufficiently complete form. 


The present book is an attempt to satisfy the above-mentioned 
demand. It should be understandable to all those acquainted with the 
mathematics courses of the higher technical institutes as well as to 
students in the senior courses of these institutes. 


A small list of the basic literature on astronautics and the classical 
two-body problem existing in Russian is presented at the end of the book. 


In writing the book the author used material from lectures on the 
mechanics of space flight intended for students of the Lomonosov University 
of Moscow, with the kind permission of D. E. Okhotsimskii, Corresponding 
Member of the USSR Academy of Science. Some of the material was used 
in the explanation of sections 3, 4, 7, 10 and 11. Moreover, comments and 
remarks of Docent I. I. Razin and Engineer B. F. Zhdaniuk, who attentively 
read the book in manuscript, contributed to an improvement in the 
explanation. In this connection, the author takes this opportunity to extend 
to them his full gratitude. 


2 April 1961 D. Pogorelov 
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INTRODUCTION 


Among the many problems of near-earth and cosmic flights, one 
problem occupies a very special place: that of studying the free motion of 
a flying vehicle in such proximity to the celestial body under investigation 
that all the forces acting on the freely flying vehicle, except the force of 
attraction produced by the celestial body, can be neglected. The peculiarity 
of this problem, in contrast to the majority of other problems of near- 
earth and cosmic flights that present great difficulties to investigators, 
is that it reduces to the well-studied two-body problem of celestial 
mechanics and presents no great difficulties as to mathematical de- 
scription andstudy. Hence, all the motions included within the scope of the 
two-body problem are usually called Keplerian motions (after J. Kepler, 
who discovered the fundamental laws of the elliptical motion of the planets, 
which are also valid for a wider class of motions), Thus, Keplerian 
motions are the simplest kinds of motions of flying vehicles in cosmic 
space. It is known that the trajectories of the more complex space flights 
can be represented, in a first approximation, as consisting of Keplerian 
trajectories. Thus, for example, the flight trajectory to any of the planets 
of the solar system can be represented as consisting of an initial Keplerian 
orbit with a focus at the center of the earth, an intermediate Keplerian 
orbit with a focus at the center of the sun, and a final Keplerian orbit with 
focus at the center of the planet to be reached. 


It is hence natural that the theory of space flights should start with the 
theory of Keplerian motions. 


In the scheme of Keplerian motions, the flying vehicle is taken as a 
point mass moving in a plane passing through the center of the celestial 
body producing an attractive force. The study of the properties of 
Keplerian motions and the establishment of fundamental dependences 
between their characteristics when analyzing the motion of flying vehicles 
in this plane is the basic content of the present hook. We start the ex- 
planation with the formulation of the requisite theorems of the mechanics 
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of a point mass, the determination of the attractive force of the celestial 
body which has a spherical mass distribution configuration, the formula- 
tion of the classical two-body problem, and the derivation of the trajectory 
equation for Keplerian motions. Then a derivation will be given of the 
relations connecting the trajectory parameters to the initial motion con- 
ditions and the magnitude and direction of the flight velocity; the angular 
range and the flight time for any current point of motion are determined. 
Furthermore, we analyze the family of trajectories of identical initial 
velocity and identical range and we determine the maximum range and the 
minimum initial velocity trajectory parameters. Finally, at the end of 
the book relations will be derived which connect angular range and flight 
time deviations with initial motion deviations. 


We shall take the flying vehicle to be a free point mass. It will hence 
be convenient to start the explanation not with an analysis and integration 
of the equations obtained by projecting the original vector equation 
describing the point motion onto the axes of some coordinate system, but 
with an application of fundamental theorems of mechanics of a free point 
mass to the kind of motion to be studied. 


1. FUNDAMENTAL THEOREMS ON THE MECHANICS 
OF A FREE POINT MASS 


The direct integration of the equations of motion of a free point mass, 
which are obtained by projecting the initial vector equation 
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onto the axes of some coordinate system, where m is the mass of the point, 
r the radius vector of the point, F the total force acting on the point, and’ 
t the time, is often very difficult. In certain particular, but important, 
cases these difficulties can be eliminated by using a transformation of the 
initial vector equation (1.1). The general dependences obtained by the 
transformation of the initial equations of motion of the point, and which 
permit the first integrals of the initial equation to be obtained in certain 
particular, but often encountered, kinds of forces are called the general 
theorems of the mechanics of a point mass. 


The problem under consideration below refers directly to sucha 
case, when it is convenient to use theorems on the kinetic moment and on 
the kinetic energy. 


1. Kinetic Moment Theorem or the Theorem of Areas, This theorem 
is obtained by the following simple transformation of the original equation 
(1.1). By multiplying both sides vectorially by r , we obtain 


a ifn ee 
Fe MG hE 
dt 


after which, by using the identity 


This dependence is called the kinetic moment theorem and is formulated 
as follows: The time derivative of the kinetic moment of a point mass 
relative to a certain center equals the moment, with respect to the same 
center, of a force under the action of which the point moves. 


This theorem is often given another form also. It is known that the 
quantity 


d> = r x dr 


is the oriented, doubled area of a triangle formed by the vectors r. dr 
and r + dr, or equivalently, the oriented, doubled area of a sector 
described by the radius vector r inthe time dt (Fig. 1). The derivative 


(1.3) 72 225 Fe = 


is called the doubled sectorial velocity. 
By expressing the left side of (1.2) in 
terms of the doubled sectorial velocity, 
we obtain 


Fig. 1 (1.4) ma = rx F, 


Since the sectorial velocity expresses the velocity at which a certain 
area is swept by a radius vector, the kinetic moment theorem is called the 
theorem of areas in later treatments. It can be formulated in the following 
manner: The product of the mass of a material point by its doubled sec- 
torial acceleration, defined with respect to a certain center, equals the 
moment, with respect to the same center, of a force under whose action 
the point moves. 


If the moment of the force acting on the point mass equals zero, with 
respect to a certain center, during the whole time the point moves, then 
the kinetic moment with respect to this center remains constant on the 
basis of (1.2): 


(1.5) rxm—= ec. 


A force with this property is called a central force; the center relative 
to which the moment equals zero is the center of the force; the motion 
performed by the point under the action of the central force is central 
motion. 


Hence, in the case of a central force the kinetic moment theorem 
yields one vector or three scalar first integrals of the differential equa- 
tions of motion of a point mass. 


According to (1.3), the vector integral (1.5) can also be represented 
in another form: 


(1.6) oF 
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which shows the constancy of the doubled sectorial velocity @ in the case 
of a central force. The vector integral (1.6) and the three corresponding 
scalar integrals are called area integrals. 


2. Theorem of the Kinetic Energy or Vis Viva Theorem. This theorem 


is also obtained simply by using a transformation of the original equation 
(1.1). 


Let us write the obvious equality 
(1.7) — dt = dr. 


Multiplying the left and right sides of (1.1) and (1.7), we obtain: 


— — 


m(v - dv) = F- dr 


Finally, let us write the desired dependence thus: 


2 
(1.8) a (me") = F.dr. 


This dependence is called the kinetic energy theorem or the vis viva 
theorem. It is formulated as follows: The differential of the kinetic energy 
of a point mass equals the element of work by a force acting on a point. 


If the force acting on the point mass has a force function U and, 
therefore, the element of work can be expressed by the total differential 
of the force function 


F. dr = dU ’ 


then (1.8) can be integrated directly and leads to the first integral of the 
equation of motion of a point 


where h is an arbitrary constant. This integral is called the energy or 
the vis viva integral. 


2.ATTRACTIVE FORCE OF A MATERIAL BODY 
HAVING SPHERICAL STRUCTURE. 
BASIC EARTH ATTRACTIVE FORCE 


According to the law of universal gravitation, also known as Newton's 
law, each two material particles are mutually attracted with a force 
proportional to their masses and inversely proportional to the square of 
the distance between them. 


Let M and m be the masses of two material particles, and r their 
mutual distance; then each particle is acted upon by a force directed 
toward the other particle which equals 


(2.1) F 


| 


where f is the universal gravitational constant. 


Often only the attractive force acting on one particle, the particle with 
mass m, say, must be examined. The effect of the particle of mass M 
on the particle of mass m causes this force, and when m = 1 


_ iM 
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It is known from mechanics that the previously mentioned force is con- 
servative, i.e., has a force function. The force function of the force (2.2) 
is called the Newtonian potential of the particle of mass M and it equals 


(2.3) Ve= i= 


where r should be considered a spherical coordinate of the particle of 
mass m = 1 in a coordinate system in which the particle of mass M 
is taken as the origin. * 


It follows from the properties of the force function that the projection 
of the force (2.2) on an arbitrary x-axis is 


(2.4) Fy. = = 


In particular, the projection of the force (2.2) on the r-axis directed from 
the particle with mass M toward the particle with unit mass equals 


(2.5) Fo == 


The physical meaning of the potential (2.3) is easily expounded by using 
(2.5). 


*From the viewpoint of the terminology used in theoretical mechanics, this 
function should be called the Newtonian work function. However, the somewhat 
inaccurate designation "potential" has been inculcated in the extensive literature 
on the theory of the Newtonian attraction, and we shall also retain it. 


Let us consider a force F' which is equal in magnitude to the force 
F with which the particle of mass M attracts the particle of unit mass 
but is opposite in direction. Evidently 


av 
, ee on il 
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Let us apply the force F' toa particle of unit mass and let us 
determine the work performed in displacing this particle from a certain 
initial position, characterized by the quantity r, to infinity. Asis 
shown in courses on theoretical mechanics, the work performed by a 
conservative force is independent of the path and is determined completely 
by the positions of the initial and final points. Hence, the desired work 
will be independent of the trajectory of the displaced particle because of 
the conservative nature of the force F' and, for simplicity of the 
calculations, the particle can be removed to infinity along a straight line 
passing through the particle of unit mass and the particle of mass M. 

In this case the work to be determined is 


(2.7) <A = fr dr = V(r) - Vio) = V(r) 
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From (2.7) it follows that the Newtonian potential at some point 
equals the work which must be expended to overcome the force of 
Newtonian attraction when removing the attracted particle of unit mass 
from that point to infinity. 


Just as the Newtonian potential of a single particle has been intro- 
duced above, it is also possible to introduce the Newtonian potential of 
n particles: 


as 
(2.8) V= ie 
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where M, is the mass of the i-th particle, r, the distance of the i-th 
particle from the particle of unit mass. 


Correspondingly, the Newtonian potential of an arbitrary body of 
mass M can be expressed in the following general form: 


(29) var f@ 
& 
M 


where r is the distance of the body particles which have mass dm 
from the particle of unit mass, and the integration is taken over the 
entire mass of the body. 


The problem of determining the integral (2.9) for a body of complex 
shape and complex mass distribution is very difficult. However, it is 
easily solved for bodies having a spherical structure. 


A body which has the shape of a sphere and for which the density at 
all points equidistant from the center is identical will be called a body of 
spherical structure. 


Let R denote the distance of an arbitrary point of a body of 
spherical structure from its center and p_ the density at this point. Then, 
evidently 


p = p(R) 


Let us describe a spherical surface of radius R from the center of 
the body A (Fig. 2). 
Furthermore, taking as axis 
the line AO connecting the 
center A to the particle 
being attracted, O, let us 
describe two conical surfaces 
ABB! and ACC' with the 
common vertex A and 

Fig. 2 vertex angles 2g and 

2p + dg), respectively. 

The surface area of the spherical band BCC'B' cut by these cones will 
equal 


27R sin go R do , 
Furthermore, let us describe an analogous spherical surface of 
radius R + dR. The mass of the part of the body enclosed between 
these two spheres and two cones will evidently equal 


dm = /— 2nR“sin gy GR dg . 


The distance r of all the particles of the part of the body thus isolated 
from the attracted particle O can be considered identical. Then 


R 
0 T 
sin d 
v= f %™ = ot f pm Rear f “2 
M * 0 0 


where Rg is the radius of the whole body under consideration. 
Setting the distance AO equal to rg, let us write 
Ou pe we ee 2 
r R* + Lo 2Rr cos Oo. 


After this, by replacing the variable of integration gy by the variable r 
which will vary between the limits LQ - R to lo + R, we obtain: 


r dr = Rr,sin gy dy 


and 
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Since the mass of a body of spherical structure is 


then 


(2.10) Ve=f 7 
"0 
The expression (2.10) shows that a body of spherical structure 


attracts an external point just as if all its mass were concentrated at the 
center. 


Since the earth, like other celestial bodies, has a very complex 
structure which has been studied only slightly, it must be replaced by 
some simpler model in many cases. In the simplest model used, the 
earth is represented as a body of spherical structure. In this case the 
Newtonian potential and the earth's attractive force can be represented 
as follows: 


TI 
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where m is the mass of the particle being attracted; r the radius 
vector of the particle drawn from the center of the earth; 7, = f{M the 
product of the gravitational constant f by the mass of the earth M. 
The numerical value of dis constant can be taken as mg = 3.986:1014 
m/sec according to the cata of many works. , $.%ee sor is/2 
Ler om and vor 
For brevity. we shall designate the potential and the force defined ' 
above as the fundamental rctential and the fundamental force of earth's 
attrection. 


The acceleration communicated to the particle being attracted by the 
fundamental force of earth's attraction is 


(2.183) 6 =- @F 


We shall call this the acceleration of the fundamental force of earth's 
attraction. 


3. CLASSICAL TWO-BODY PROBLEM. EQUATION OF 
THE TRAJECTORY OF A POINT MOVING UNDER 
THE EFFECT OF THE FUNDAMENTAL FORCE OF 

EARTH'S ATTRACTION 


Let us assume that there are two bodies which are mutually attracted 
according to Newton's law. Let us also assume that they have a spherical 
structure. Then on the basis of the exposition in the preceding section, 
the attraction between these bodies can he considered as coincident with 
the attraction of two points having the masses of these bodies and located 
at the centers of mass of these bodies. Therefore, the problem of 
studying the motion of two bodies under the effect of their mutual 
attraction, which is called the two-body problem in mechanics, reduces, 
for bodies having spherical structure, to the considerably simpler 
problem of studying the motion of two points under the effect of their 
mutual attraction. The simple prohlem involving spherical bodies is 


usually called the classical two-body problem, to distinguish it from the 
general two-body problem. In those works where only this problem is 
analyzed and, therefore, there is no necessity to distinguish it from more 
general problems, it is usually called simply the two-body problem. As 
is well known, it was first solved by I. Newton by a geometric method. 

D. Bernoulli gave its analytical solution, and finally, it was investigated 
in detail by L. Euler. Up to now, changes in its solution have involved 
mainly the choice of variables in terms of which it will be expressed. 

A detailed explanation of this problem is usually given in all courses on 
celestial mechanics and theoretical astronomy. 


The motion of a projectile or any apparatus launched from the earth's 
surface into a vacuum space can, in many cases, be reduced to the scheme 
included within the scope of the classical two-body problem. In particular, 
the problem of studying the motion of a projectile or any apparatus ina 
vacuum space near the earth can be represented as a classical two-body 
problem in which one body has a spherical structure and mass M, equal 
to the mass of the earth, and the other can be considered as a material 
point having mass m equal to the mass of the projectile or apparatus. 
Let us examine this particular case of the classical two-body problem 
for which mass M is so much larger than mass m that it is possible to 
neglect the effect of the mass m on the mass M. 


Let A denote the center of mass of the earth and S that of the 
point mass whose motion we shall analyze. Because of the foregoing, 
the motion of the point S can be considered as taking place under the 
effect only of the fundamental force of the earth's attraction Fy. directed 
toward the point A which will be considered fixed. The usual means 
could be used to write the equation of motion of the point S; we take as 
the original equation 


re 


dF 5 
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and by successive transformation and projection on the appropriate axes 
we obtain all the dependences of interest. However. this problem can be 
solved more briefly and simply by using the theorems of the kinetic 
moment and the kinetic energy. 


Thus, by using the kinetic moment theorem it is possible to obtain a 


vector integral of the equation of motion of the point S in the form of 
(1.5) 
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where r is the radius vector of the point S drawnfrom A and ¢c is 
a certain constant vector. 


This integral shows that the point S can perform motion only in a 
plane passing through the point A and perpendicular to the constant 
vector c. It is evident that because of this motion it is convenient to 
exumine the point S in an Axyz coordinate system whose origin is at 
the center of the earth, the Axy plane being superposed on the plane in 
which the motion of the point S occurs. In the coordinate system thus 
selected, the study of the motion of the point S under the effect of the 
fundamental force of earth's attraction reduces to the study of its motion 
inthe Axy plane. The vector integral written above thereby yields the 
simplest solution of a second-order differential equation determining the 
change in the coordinate z of the point under consideration in the form 


dz) _ 4g. = 
at 0; z(t) 0 

Moreover, since the vector integral is equivalent to three scalars, then 

as will be shown below, it will also yield a first integral of the equation 

of motion of the point S inthe Axy plane. 


Because of the spherical symmetry of the field of fundamental force 
of the earth's gravitation, it is convenient to examine the motion of the 
point S inthe Axy plane in the polar coordinates r and 7, related to 
the Ax- and Ay-axes as shown in Fig. 3. In the figure, S is the 
current position of the point under consideration and QO, is a certain 
initial position. In our case it is convenient to determine the velocity 
vector v_ of the moving point S by using the absolute value of this 
vector v and the angle of its inclination § to a plane perpendicular to 
the radius vector of this point. It is customary here to take the upward 
direction from the plane perpendicular to the radius vector as the positive 
direction of measuring § and the downward direction (toward the center 
of the earth) as the negative direction. The angle @ thus defined is 
called the angle of deviation. In conformity with the notation introduced, 
the parameters defining the position and velocity of the point S will be 
denoted thus: for the current instant, r, 7, v, 9, and for the selected 
initial instant r,, 14, V4; 04 
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We must obtain the scalar integral of the equations describing the 
motion of the point under consideration inthe Axy plane from the vector 
integral written above. Hence, we decompose the radius vector of the 
point S and its velocity vector into components along the polar 
coordinate axes: 

is +0 70 , 


= —~ = med @) 
(3.1) y = v,r Vi, : r= rr 


where r9 and 7? are unit vectors of the polar coordinate axes drawn 
through the point S (Fig. 3). 


Fig. 3 


Using (3.1), let us transform the vector integral mentioned; we 
obtain: 


rxmy * m [rr x (Vpre + v,7 | 
= mrv, (x? x 79) =c. 
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Taking into account that 
(3.2) V2.3 Vv. = gol 


and that the vector product of the unit vectors r9 and 79 is a unit vector 
parallel to the constant vector c and 


Cc es 
— = 0 = const, 
ni 


the desired scalar integral can be represented on the basis of (1.6) in the 
form 


(3.3) r2 21 = 0 = const. 


Equation (3.3) shows that the magnitude of the doubled sectorial 
velocity o of the point S with respect to the center of the earth remains 
constant throughout its motion. 


However, knowledge of the integral (3.3) is still insufficient to 
determine the trajectory of the point S. The trajectory equation can be 
written by using, besides this integral, the kinetic energy theorem in the 
form (1.8): 


i 2 _F .- 
(3.4) 5 d(v*) = — dr. 


In the case we consider, on the basis of (2.12), 


F “) 
jee ge 0 
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re 


(3.5) 


Substituting (3.5) into (3.4), we rewrite it thus: 


1 aqw*) _ _ 0 =o ar 
2 dt 2) dt 


Using (8.1) for the velocity vector and relations (3.2), let us express 
the equality of interest in terms of scalar quantities only: 


13 


2 dt re at 
This equality relates the three variables r, v and t. In order to write 
the trajectory equation in the polar coordinates r and 7), it is necessary 
to obtain an equation relating only these two variables. Such an equation 
can be obtained from (3.6) if a transformation is made from the variable 
t to the variable » and then v is expressed in terms of r and 7. 


Let us first transform from the independent variable t to the 
independent variable 7 in (3.6). To do this, we rewrite it as follows: 


This equality is equivalent to the two equalities: 


ay 
(3.7) = 0 
and 
aig, ae 
2 dn ye an 


The first equality permits us to obtain the trajectory equation directly 
for the simplest particular case of point motion, the case of rectilinear 
motion 

(3.9) n= Io = const . 

The second equality permits the elimination of the variable t and, later, 
of the variable v and permits us to obtain the trajectory equation for the 


general case of point motion. 


In order to express the quantity v intermsof 7 and r, we use 
the expression 


2 2 
2 _ (dr 2 (dy 
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which can easily be obtained from (3.1) and (3.2). Representing (3.10) in 
the form 


2 2 
9 _ {dy dr 2 
(3.11) v% = (22) () + 2 


and using the integral (3.3), we eliminate the variable t: 


o2 |/ar\? D) 
(3.12) v -= (2) + Yr 


Then, substituting (3.12) into (3.8), we obtain the desired trajectory equa- 
tion. However, before this it is convenient to replace r by the new 
variable 


(3.13) u = 


re Le 


which can be used to write the differential equation of the trajectory in its 
simplest form. Transforming to the variable u in (3.8) and (3.12) by 
using the expressions 


1 dr 1 du 
(3.14) re= 7 and ant “72 an ® 
we obtain 
3.15) 2 aw) _du,. v2 = o@ |()" 4 2 
(3.15) 2 dn dn 0’ dn 


Now substituting the expression for v- in the first equality, the desired 
equation can be written thus: 


T 
(3.16) au pay ee. 
dn dn og” dn 


The equality (3.16) permits the trajectory equation to be obtained for 
the particular case of point motion along a circle 


(3.17) —=0 or ~—=0 


or 
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(3.18) Yr = Lg = const. 


For the general case of point motion, on the basis of (3.16) we can 
represent the desired differential equation of the trajectory as follows: 


2 
(3.19) a + u = Ug 
dn 
where 
™ 
(3.20) U, = —= = const. 
0 52 


The general solution of the linear differential equation (3.19) with 
constant coefficients can be written as:* 


(8.21) u = Up - U,cos7 - Upsin 7 


1 


where U, and U,. are integration constants and U, is the particular 
solution of the inhomogeneous equation (3.19); or it can be reduced to the 
form: 


(3.22) u = U, E - ecos(n - ny)] 


where the constants e and nn. can be expressed in terms of the constants 
Ug, Uy, Us as follows: 


2 2 
U7 + U 
(3.283) e? =—-—_—4 
Ue 
0 
U U 
1 2 
(3.24) cosnp = | sin), = . : 
U, + 
1+ “a se 


(The geometric meaning of 7 F will be clarified below.) 


*For example, see V. I. Smirnov: Kurs vysshei matematiki (Course in Higher 
Mathematics) vol. 2, Gostekhizdat, 1951, sections 27. 28. 
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The relations (3.23) and (3.24) can easily be obtained by comparing 
the expressions (3.21) and (3.22). 


Transforming from the auxiliary variable u to the variable r in 
(3.22) and introducing the new constant 


(3.25) p= 
Up 


let us write the trajectory equation in the final form 


ee ee 
(3.26) r 1-e€ cos(n = Ta) ’ 


Comparing the obtained equation with the equation 


p 
3.27 P= 
( ) 1 - e cos B 
which is the equation of a conic section in polar coordinates,* we 
arrive at the conclusion that the trajectory of a point moving under the 
effect of the fundamental force of earth's attraction is a conic section. 


Let us recall that in 
writing (3.27), one of the 
foci of the conic section is 
taken as the pole of the 
polar coordinate system 

ELLIPSE and the focal axis of the 
conic section is taken as 
the polar axis. The line 


MYPERBOLA 


POLAR AXIS 


DIRECTRIX 


focus is chosen as the 
positive direction of the 


from the focus chosen as the 
pole drawn through the other 


polar axis for the ellipse, 
Fig. 4 and the direction from the 
focus chosen as the pole to 


infinity on the branch of the curve under consideration is taken for the 
parabola and hyperbola (Fig. 4). 


*For example, see I. I. Privalov: Analiticheskaya geometriya (Analytic Geometry), 


Gostekhizdat, 1945, Chapter 4, section 4. 
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Let us clarify the geometric meaning of the parameters in the tra- 
jectory equation. As is known, a conic section can be defined as the 
geometric locus of points, the ratio of whose distances from a given point 
(focus) to the distance from a given line (directrix) is a constant value, 
equal to e. Hence (see Fig. 4): 


‘ | FM, 
or the ellipse e = aN. <1; 
11 
FM, 
for the parabola e = —— = 1 ; 
M,N 
22 
FM, 
for the hyperbola e = —— > 1 
M3N3 


Therefore, the parameter e determines the shape of the curve. It is 
called the eccentricity of the conic section. 


As is seen from (3.27), the quantity p equals the length of the radius 
vector of a curve drawn from the focus and perpendicular to the focal axis. 
Hence, the parameter p defines the linear dimension of the curve. It is 
called the focal parameter of a conic section. 


As is seen from (3.26) and (3.27), the quantity Ny equals the angle 
formed by the focal axis of the conic section and the polar axis of the 
coordinate system selected for the writing of equation (3.26). Hence, the 
parameter 7, defines the position of the focal axis of the conic section 
relative to the selected polar axis. 


Summarizing, we can make the following deduction: The trajectories 
of points moving under the effect of the fundamental force of the earth's 
attraction can, depending on the initial conditions, be radial, circular, 
elliptical, hyperbolic or parabolic. For brevity, all these trajectories are 
called Keplerian trajectories and the motions pertormed by the point mass 
along these trajectories are Keplerian motions. 
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4. RELATION OF THE FUNDAMENTAL PARAMETERS 


OF KEPLERIAN TRAJECTORIES TO THE INITIAL 
MOTION CONDITIONS 


In order to determine the influence of the initial conditions of the 
point motion on the shape, magnitude and location of the trajectory, it is 
sufficient to write the relations connecting the fundamental trajectory 
parameters e, p and 1), with the point-motion parameters corresponding 
to a certain instant taken as the initial instant. Since any arbitrary instant 
of motion can be selected as the initial instant, we will not make it definite 
here but will write the desired relations in a form valid for any time of the 
point motion under examination. 


Let us introduce the nondimensional point-motion parameter 


= 2 
44.1) k= = reas 
"0 "O 


which is the ratio of the kinetic energy of the point at any instant of the 
motion to the work which must be expended to transport this point to 
infinity from the position it occupies at the instant under consideration. 
This physical interpretation of the parameter k follows from (2.11) for 
the earth's Newtonian potential and from (2.7) which is valid for any 
Newtonian potential. The convenience of taking this parameter into con- 
sideration is explained by the fact that the trajectory shape is completely 
defined by its magnitude alone. We arrive at this conclusion somewhat 
later, but it can naturally be expected by starting from the physical 
meaning of the parameter. 


The motion under consideration can be defined completely if the 
coordinates of the moving point and its velocity vector are assigned at an 
arbitrary instant, say, by using the quantities 


(4.2) r, 7, v, 8 


introduced in the preceding section (Fig. 3). 
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By expressing the trajectory parameters e, p and Ny in terms of 
the quantities in (4.2), we obtain the desired relations. To solve the last 
problem, let us first express the integration constants 


(4.3) Ug, Uy, Us, 


in terms of the quantities in (4.2). Returning to (3.21), we write 


. du U 
(4.4) u = Uy - U jcos7 - Uosin7 ; an U,sinn - U,cos7 . 
Moreover, it follows from (3.20) that 
(4.5) UU, = 5 


Making use of (3.3), (3.13) and (3.14), we write: 


= ye Gl 2 ab 
e dat’ 4 y 
(4.6) dr 
dul _ 2dr _ 1 at 
dn dn r2 an 
dt 
It follows from Fig. 3 that 
(4.7) Vac v sin? ; v, = vcos 6 
and from (3.2) that 
dr ; dj _ vos @ 
. ~~ = 6 : oOo 
(4.8) at v sin + ot = 


Substituting (4.8) into (4.6), we obtain: 


Finally, substituting (4.9) into (4.5) and (4.4), let us write the following 
relations to determine the constants of (4.3): 
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BO: og Ol 


U = = 
: y? r@v2cos20 
1 
4, j = - -— 
(4.10) U,cos n + U,sin n Uy 
tan 6 
U,sin n - U,cos n a regs 


from which, by taking (4.1) into account, we obtain 


ee: eee ae 
0 = = 
of oy-e6s-6 2kr one 
1 tan é@ 
U, = (u, - 7 cos] - — sin? 
(4.11) oe a . : tan 6 
9 = 0 = sinn coOs7) 
2 2 
24 72 = il tan“0 
Uy US (v, ) - r 


Using (3.25), (3.23), (3.24) and (4.1) after this, the desired relations can 
be written in the following general form: 


rq) 
bo 
1} 
Gc 
= bo 
+ 
e 
bo bo 
“ee 
( 
ie 
ed Sw 
bo 
“oe 
et 
a) 
E 
D 
NS 
to 


(4.13) 4 e 26)? + (2k sin @ cos 9)? 


(1 - 2k cos 


(1 - 2k é65-0)- + 4k? (1 - cos“6 )cos26 


II 


1 - 4k cos“é + 4k2cos6 


from which 


2 


(4.14) e? = 1 - 4k(1 - kjcos“6 


or 
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2 
(4.15) e* = (1 - 2k)“cos6 + sin 0 


Uo (Ugr - 1)sinn + tan 6 cos 7 


Gas) hee = a (Ujr - Icosn- tan é sing 


After having written the desired expressions, it is possible to write 
the dependence of the shape, magnitude and location of the trajectory on 
the initial conditions of the motion of the point under consideration. 


Let us determine the condition guaranteeing point-motion along the 
rectilinear trajectory for which we obtained (3.7) or (3.9) in the previous 
section. It follows from (3.7) and (4.8) that motion along a rectilinear 
trajectory is possible only when there is compliance with the equality 


(4.17) 06 = 


bo 13 


even at any initial instant of the motion. 


Now, let us determine the condition guaranteeing point-motion along 
a circle. The trajectory equation for this case was written in the 
preceding section in the form of (3.17) or (3.18). It follows from (3.17) 
and (4.9) that motion along a circular trajectory is possible only when the 
equality 


(4.18) 6 =0 


is always satisfied, even at any initial instant. However, compliance with 
condition (4.18) at some initial instant, which is the necessary condition, 
is not sufficient for circular motion, since it does not guarantee zero 
eccentricity. The additional condition, which must be satisfied at the 
selected initial instant for the trajectory to be circular, is also easily 
obtained from (4.15), with (4.18) taken into account: 


Ee = (1 - 2k)"cos“6 + sin’ 6 = (1 - ok) = 0, 


from which 


2 
rv 1 
(4.19) kk “= On, = ° ‘ 


—) 
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Evidently this condition must be satisfied at any instant for motion 


along a circle, since any instant can be taken as the initial instant of the 
motion. 


In all the remaining cases, the trajectory shape is determined com- 
pletely by the magnitude of the parameter k and, as is easily seen, from 
(4.15): 


Always e < 1 for 0 < k < 1, i.e., the trajectory is an ellipse; 


1 for k 


Always e 1, i.e., the trajectory is a parabola; 


Always e > 1 for k > 1, i.e., the trajectory is a hyperbola. 


This result corresponds completely to the physical deductions which can 
be made by starting only from the physical meaning of the parameter k. 


In conclusion, let us recall that a very convenient parameter for 
elliptical trajectories is also the semi-major axis a which is expressed 
as follows in terms of the parameters p and e: 


Pp 
(4.20) a= {= ee 


Substituting (4.12) and (4.13) into (4.20), we obtain 


, 
(4.21) a= oa, soy * 

If it is agreed to consider the semi-major axis for the hyperbola to 
be negative, then (4.20) and (4.21) will be identically applicable for both 
the ellipse and the hyperbola and we shall use them repeatedly later. 


In addition to the fundamental parameters we have examined in this 
section, we will be interested in many other parameters of Keplerian 
trajectories. These parameters and their relation to the initial motion 
conditions will also be examined in succeeding sections. 
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5. MAGNITUDE OF THE FLIGHT VELOCITY 


The magnitude of the flight velocity at any arbitrary point of the 
trajectory can easily be expressed in terms of the initial motion con- 
ditions by using the kinetic energy theorem as written in the form (3.6). 
Rewriting (3.6) thus: 

20 
(5.1)  d(v2) = -—~ adr 
Yr 


we obtain after integrating 
(5.2) v’ = v2 4+ 27(2 z +) 


where v and r are the magnitudes of the velocity vector and the radius 
vector of a point at any arbitrary instant of the motion, and ‘eT and Ty 
are the same quantities at the initial instant chosen. 


Later, the relation (5.2) will often be used, reduced to the form: 


2 2 - 
(5.3) vo=v, + 2h, r, 
where 

— hy 
(5.4) hep hp try - es 


or to the form 


(5.5) k= 


where according to (4.1) 


2 
(6.6) k=—;k, = 
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As follows from the derivation, the above formulas are valid for all 
kinds of trajectories. 


If the trajectory of the point under consideration terminates on the 
surface of the earth (which we shall consider to be a sphere of radius R 
consistent with the assumptions of section 2 - see Fig. 3, for example) 
then the magnitude of the velocity at the point of intersection of the 
trajectory and the surface of the earth is usually of interest. Let us call 
this point the incident point and denote it by Op. In what follows, we will 
ascribe the subscript R_ to all quantities at the incident point. According 


to (5.3), the magnitude of the velocity at the incident point of interest will 
equal: 


2 2 = 0 
(5.7) vp =v, + 2h 
where 
(5.8) hy =; hy, = 


and correspondingly: 


k +h 
1R 
where 
Rv, 
(5.10) kR = 2m 


6. FLIGHT ALTITUDE 


Let us designate the altitude of a point moving over the surface of the 
earth h, defined by the relation 
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(6.1) h=r-R 
as the flight altitude. 

For a point moving along a radial trajectory the flight altitude 
achieves its maximum value at the instant when the point, departing from 
the surface of the earth and gradually losing speed, stops completely. 


Letting ryax denote the maximum distance from the center of the earth 
for this case, we obtain from the expression (5.2) 


(6.2) xr -_— 


and correspondingly 


(6.3) hnax = Tmax ~ FB 
As follows from (6.2), for ky = 1 the point recedes from earth toward 
infinity: 

Tmax =~ : 


As a point moves along a circular trajectory, its altitude and speed of 
flight remain constant and equal to a certain initial altitude and initial 
velocity. As a point moves along an elliptical trajectory, its flight 
altitude varies between a certain minimum value Bnnin? Called the 
perigee, and a certain maximum value hax, called the apogee, if the 
trajectory is located completely outside the limits ot the earth (see Fig. 4). 
And in conformity with the trajectory equation (3.26): 


Pp 
hmin = min ~ R$ Tin = 1lt+e 


(6.4) 


Dmax = Tmax - Ri Tmax Lee 


If the trajectory intersects the earth's surface as is shown in Fig. 3, 
then the flight altitude varies between zero and a maximum value which 
will be denoted by H in this case and will be defined exactly as above: 


(6.5) H = - R 


In this case the quantity H is often called the trajectory altitude. 


When a point moves along a parabolic or hyperbolic trajectory, the 
flight altitude changes from a certain minimum value hyyin to the 
infinitely large h,,4, = ~ if the trajectories are located completely out- 
side the limits of the earth (see Fig. 4). And in conformity with (3.26) the 
value of h,,j, is also determined by means of formula (6.4) in this case. 


If the trajectories intersect the surface of the earth, then the flight 
altitude changes from zero to infinity in the case of a point moving from 
earth, and from infinity to zero in the case of a point moving to earth from 
an infinitely distant point. 


t. ANGULAR FLIGHT RANGE 


One of the convenient characteristics of point motion near earth is 
the flight range measured along the earth's surface from the radius vector 
of some one point of the trajectory O, to the radius-vector of some other 
point Oo. Letting ®}5 denote the angle formed by these radius vectors 
drawn from the center of the earth, the flight range along the earth can be 
characterized completely by this angle. The linear range measured along 
the surface of a spherical earth will here equal 


Gi), Dig = hoy 


In the following, in contrast to the linear range Dio: the angular 
range Do will be called the angular flight range. 


The determination of the angular range % between two trajectory 
points is not difficult if the coordinates of these points are given. The 
problem is not related just to the trajectories or the point motion under 
consideration. Naturally, the formulation of such a purely geometric 
problem as the location of two points whose coordinates are given is not 
of interest in the study of the kinds of motion under consideration here. 


Another formulation of the problem of determining the angular range 


between two points of the trajectories being examined is more expedient 
for what follows; it is the following. 
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The parameters (4.2) 


(7.2) r , V4, @ 


Nye Vy 4 


’ 
which completely define the trajectory, are given at one of the trajectory 
points, say Oj. Itis required to determine the angular range Dio 
between the point O, anda point O» of the trajectory, which has only 
one coordinate lo given. By giving the point O»5 just the one coordinate 
Io, we have posed the problem of determining the angular range hetween 
the point O,, at which the parameters (7.2) have been given, and the point 
of trajectory intersection with a sphere of radius lo described from the 
center of the earth. Hence, this problem can also be represented as a 
problem to determine the angular range of a point of incidence on a sphere 
of radius r,. by means of prescribed initial conditions (7.2). The 
determination of the angular range here reduces to the expression of the 
angle ®,. in terms of the initial conditions (7.2). 


Evidently, in the case of point motion along a radial trajectory the 
angular range between any of its points is zero. The formulated problem 
has no solution for point motion along a circular trajectory: for Ty # lo 
it is missing and for ty =, it seems indeterminate. 

In the case of an elliptical, parabolic or hyperbolic trajectory, the 
problem of determining the angular range can be solved as follows. 


In order to simplify the subsequent calculations, let us introduce a 
new polar axis passing through the center of the earth and the prescribed 
point O,, and let f denote the polar angle of the moving point, measured 
from the new polar axis which is connected to the prescribed point O; 
(see Fig. 3). Then, evidently 


(7.3) f=n - 7 
and the new polar angles of the points O01, and Oo will equal: 
(7.4) f, =O, f, = = © 

Since the trajectory equation and all the preceding formulas have 
been written for any arbitrary position of the polar axis passing through 
the center of the earth, they remain valid even for the special direction 
selected above of the polar axis connected to the point O,. Hence. the 


trajectory equation in the new polar coordinates can, as before, be 
described in the form of (3.21) 
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7.5 = - = 
(7.5) u Uy U,cos f U,sin f 


and the constants U,, U1; U, which enter here, can be defined as 
before by means of formulas of the form of (4.11): 


1 
0 2kr cos26 


ae 
i 


(7.6) 4 U, (u, - *) cos p - OS ging 


= a5 oe tan 0 
U, = (vo - *) sin f + : cos f 
and by the parameters 


(7.7) Ta f, = 0; Vy o, 


prescribed at the point O,. After this, the problem of determining the 
angular range reduces to the determination of the point of intersection of 
the trajectory (7.5) with a sphere whose equation has the form 


1 
(7.8) u=uUu, =>, = const . 


It is evident that to determine the desired point of intersection, it is 
sufficient to substitute (7.8) into (7.5) and, by solving it with respect to 
{, to obtain the unknown coordinate of the point Oo: 


(7.9) f=f, = ® 
which equals the desired angular range, whereupon the formulated 


problem is solved. 


Thus, by substituting (7.8) and (7.9) into (7.5), we obtain the initial 
equation for the determination of the angular range in the form 


= - - in ®& 
(7.10) uy, Uy U,cos ® U,sin 


where, in order to simplify the writing, it has been assumed that 


(7.11) B= ¢. 
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Transforming to the half-angle, let us rewrite equation (7.10) thus: 


ae ee ee — ae 
(Up - Up) - U; (0s 9 sin | 2U,sin 9 cos 5 : 


& 
Dividing both sides of the equation by cos” re as is possible for 


} 
(7.12) cos” a #0 


we obtain 


or 


2: 2 py fan a U.) = 0 
+ U,)tan - an = + (Up - Ug - p= 0- 


(7.13) (Uy - u ; , 


By substituting the prescribed values (7.7) into (7.6) we obtain the 
following simple expressions for the constants Uo: Uy and Uo: 


tan 6 

1 1 1 

(7.14) U, = ———-; U, =U, -—; U, = 
0 2k jr 4cos?0, 1 0 2 r 


Using (7.14), let us rewrite (7.13) thus: 


tan @ 
& ® 
(7.15) (2v, ee =, ) tan? 5-2 Stan 8 oe = 
ry lo ry 2 Ty lo 
or 
2 2 2 
(7.16) X,tan , .2 tan a, tan a hy = 0 
where 
(7.17) a ar eee r xX, = d 2-h 
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Solving (7.16), we obtain the desired expression for the angular range 
in the following form: 


6 + /tan2 h 
tan 0, tan“@, + X,h, 


Xy 


(7.18) tan 


bo |e 
| 


Let us examine the solution just obtained. It shows that in the general 
case, depending on the magnitude of the radius r, of the sphere, we can 
expect the presence either of two points of intersection of the trajectory 
with the sphere or of one tangent point or of no common points. Let us 
first determine the radius of that sphere which is tangent to the trajectory. 
To do this, it is sufficient to equate the discriminant of our quadratic 
equation to zero. Hence, to shorten the computations, we return to the 
equation written in the form (7.13) and obtain the desired relation as 


f _ yt e+ ye 


(7.19) (u 1 9 


- U 


2 0 


from which it follows that 


2 2 
uy = Uy + U, + Uz 
or 
whi 
Uo 
1 


or in conformity with (4.12) and (4.13) 


(7.20) rp = 


By comparing (7.20) with (6.4), the deduction can be made that only 
those spheres which pass through the vertex can be tangent to the 
trajectory. 
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In particular, an elliptical trajectory has two vertices corresponding 


to the values r,i, 29d YVpyax: therefore, two spheres with the radii 
(7.21) ro => Tmin * [Tae and ry Ty te 


can be tangent to it. 


The parabolic and hyperbolic trajectories have just one vertex 


corresponding to the value ee hence, just one sphere of radius 


(7.22) Yo = Ymin = 


can be tangent to each of these trajectories. 


After the spheres which are tangent to the trajectories are determined, 
the rest are easily determined. Thus, it is evident that the elliptical 
trajectory will be intersected in two points by all spheres whose radii lie 
in the interval 


(7.23) Ds < Lo < Ene 


and will not be intersected byany spheres whose radii satisfy the following 
inequalities: 


Yr <r ., andr or 
2 min 2 max 


The parabolic and hyperbolic trajectories will be intersected at two 
points by all spheres whose radius is greater than the radius of the 
tangent sphere 


(7.25) V9 7 rmin 


and will not be intersected byany spheres whose radius is less than the 
radius of the tangent sphere 


(7.26) a) < lmin 


Let us now determine the angular range of the tangent points. On 
the basis of (7.18), the equality 
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(7.27) tan“6, + X,h, = 0 


can be written and the simple expression 


(7.28) xy = 7 = 


follows for \ 


1 Using (7.27) and (7.28) we can reduce (7.18) to the simple 
form 


°p 
(7.29) tan > = - 


1 Yo 
min max 


for the elliptical trajectory and 


ie T, ~ min 
(7.31) bh, = —— 


min 
for the parabolic and hyperbolic trajectories. 


Moreover, since the tangent points are the vertices of our trajectories 
and lie on the focal axis, then (4.16), which determines the position of the 
focal axis relative to the arbitrarily selected polar axis, can be used to 
determine their angular range. The values of U, and U, are determined 
by (7.14) for the selected polar axis passing through the given point O,- 
Substituting them into (4.16), we obtain the following expression for the 
angular range between the point O, and points lying on that half of the 
focal axis on which the second focus of the elliptical trajectory and the 
infinite point of the parabolic and hyperbolic trajectories lie: 


2k,tan 4 
(7.32) tan®, = ——, 
- 1 - 2k, + tan’ 6, 
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Fig. 5 


To determine the angular range 
of the points of intersection of the 
elliptical trajectory with the sphere, 
let us use Fig. 5. Let us first 
examine the case when 


(7.33) tan Gy > 0 and r, >To, 


that is, when the point O, is onthe 
ascending branch of the trajectory 
above the prescribed sphere. One 
of the two existing points of inter- 
section in this case can be considered 
a point of incidence on the sphere 
since it lies on the path of the 
motion of the point under considera- 
tion and it can actually be reached 
for motion from the prescribed 
point O,. Let us denote this point 
of intersection by Oo. It is possible 
to consider another point as actually 
achievable for the motion of the 
considered point from the pre- 
scribed point O;. Let it be denoted 


by O. Correspondingly, let © denote the angular range to the point O,, 
and @' to the point Oo. Then, by taking account of the positiveness of | 
tan 9, for the case considered and that 


(7.34) @' < 0, @ > 0, 


we obtain 

® 
(7.35) tan —- = 

2 

ot 
(7.36) tan 2 = 


For the case when 


(7.37) tan 6 


1° 9 and r 


we will have one of the two points of intersection more distant from the 
prescribed point O,. Let us denote it by Oo und the range to it by ©. 
Let O, denote the other, less distant point and ®' the range to it. Then, 
taking into account, as in the preceding case, the positiveness of tan 6 
and that for the new case 


1 
(7.38) @ > & > 0 
we obtain the same formulas (7.35) and (7.36) to determine © and ®'. 


For the case when 


(7.39) tan @, < 0 and Yr, 7 Yo, 


that is, when the point O, is on the descending branch of the trajectory 
(in Fig. 5 it is denoted by O; in contrast to the previous position), we 
will have the same points of intersection O, and Oo as in the first 
case. Denoting the angular range to the point Og by ® and to the point 
O5 by © and taking account of the negativeness of tan 6, and that for 
the case under consideration 


(7.40) $' <0, S$>0, ([#l > [4 
we obtain the same formulas (7.35) and (7.36) to determine © and ®'. 


Finally, for the case when 


(7.41) tan@, < 0 and r, < ry 


we will have the points of intersection 0 and O5 . Letting & denote 
the angular range to the point O, and @' the angular range to the point 


O5 and taking account of the negativeness of tan 4, and that for the 
case considered 


(7.42) <0, &' <0, [el > [4 
we obtain the same formulas (7.35) and (7.36) to determine ® and 4". 


The following deduction can be made as a result of the analysis 
carried out. 


30 


For each pair of points of intersection of the elliptical trajectory with 
the given sphere which are symmetric relative to the focal axis, the angular 
range determined by (7.35) corresponds to a point reached by moving away 
from the axis in the motion direction, and the angular range determined by 
(7.36) corresponds to a second point reached by moving away from the 
focal axis in a direction opposite to the motion. 


For the parabolic and hyperbolic trajectories cases, perfectly 
analogous considerations 
apply and an analogous 
conclusion is obtained. 

~ However, since the vertex 
PP, analogous to the vertex O, 
ue of the elliptic trajectory 
recedes to infinity for these 
trajectories and the motion 
from one point of inter- 
section with the sphere to 
the other, symmetric 
point of intersection can only 
be realized via the vertex 
for" oo (Fig. 6); then the motion 
= over these trajectories 
comparable to the motion 
over the elliptical trajectories 
considered above seems to 
be reversed with respect to 
Fig. 6 the points in which we are 
interested. The motion is 
directed from the point O, to the point O, on the elliptical trajectories 
and from the point O, to O; on the parabolic and elliptical trajectories. 

Hence, for each pair of intersection points of the parabolic or hyperbolic 

trajectories with the given sphere which are symmetric relative to the 

focal axis, the angular range determined by (7.36) corresponds to the point 
receding from the axis in the direction of motion, and the angular range 
determined by (7.35) corresponds to the second point receding from the 
focal axis in the direction opposite to the motion. 


In the sequel, we shall be interested in those points of intersection of 
the trajectory with the sphere which can be achieved by the moving point 
as it moves further from the point O, . For brevity, let these points be 
designated points of intersection lying in the direction of motion or, 
simply, forward-lying points of intersection and let the rest be designated 
backward-lying points of intersection. 
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For elliptical trajectories over which the moving point can be dis- 
placed and can complete full revolutions without encountering any 
obstacle in its path, both points of intersection with the sphere are 
forward-lying. Here, inthe case r, > r, the point O, will always be 
first in the path of the motion, followed by the point Oo (see Fig. 5). 
Hence, it is necessary to use (7.35) when determining the angular range 
to the first point in this case, and, therefore, (7.36) when determining 
the angular range to the second point. In the case ry < Yo, the point 
Os will always be first in the motion path and Ov will be second. In 
this case it is necessary to use (7.36) when determining the angular 


range to the first point and (7.35) when determining the angular range to 
the second point. 


For parabolic and hyperbolic trajectories in the case r, > ro, both 
points of intersection are forward-lying for 9, < 0 and backward-lying 
for 9, > 0. Hence, the point O, will always be the first in the motion 
path among the forward-lying points and the point O5 the second (see 
Fig. 6). Consequently, (7.35) must be used in determining the angular 
range to the first point and (7.36) to the second point, just as for the 
elliptical trajectories. We shall not investigate the determination of the 
angular range to the backward-lying points. In the case ry < Yo just 
one point, namely O5 , will always be forward-lying and the other, Oo, 
will be backward-lying, and the angular range to the forward-lying point 
will be determined by (7.36). 


Hence, in the case r lo for all trajectories, the angular range 
to the first of the forward: Pare Stuts of intersection is determined by 
(7.35) and inthe case r,; < r, the angular range to the first of the 
forward-lying points is determined by (7.36). 


8. DIRECTION OF THE FLIGHT VELOCITY 


The direction of the velocity at any point of the trajectory determined 
by the angle & can be related to the initial conditions (7.2) prescribed at 
some point O,, by means of the following simple relation: 


= = ) 
(8.1) o r4V{cos a, rv cos 


which results from (4.9). 
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In those cases when it is desirable to have a relation between the 
angles 04 and @ in terms of the tangents rather than in terms of the 
cosines, the desired dependence can be obtained as follows. Using (4.12), 
let us write 


p = 2k4r,cos76, = 2kr cos“é 
or 


(8.2) a = kr 


2 
cos @ cos 0, 


It is here evident that the cases of the radial and circular trajectories for 
which, respectively, 


(8.3) cos*@=0 and cos*9=1 
are eliminated from the considerations. 
Substituting (5.5) into (8.2), we obtain 


KT, k, + hy : 


cos“6 1 + hy cos“U, 
Making the substitution according to (5.4) 


(8.4) 1+ hp == 


we can write the last equality as 


k k xr hir- 

Pt, pee 

a. 2 
cos 9 rcos a, rjcos’9, 
or 

2 pa 
"41 1 - 

2° 27 . 

r cos @ cos eT k,cos a, 
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or 


tan*@ = tan“6, + 1- 
. 1 ee 


2,2 . . 
Replacing rr in the right side by the expression resulting from 
(8.4), we obtain 


2 = 
r h 
= Dae ni 2 - 2 1 
9 tan’@ = tan®@, + 1 - (1 + hy) + 
r kycos ¢, 


k 68-6 


2 7 1 — 
tan"@, + h -2-h, 
1 1 


Using (7.17) we can write 


2 
r 
1, 2 2 - 
— 9 = 
(8.5) 2 tan tan*6, + X,h, 


From (7.18) it follows that 


/tan® ee 2 
(8.6) £ ./tan“@, + Xzh, = X,tan> - tan 4, 


Hence, as is seen from (7.35) and (7.36), in order to obtain the positive 
root, it is necessary to substitute the value of ® corresponding to the 
points Op, and Oo in the right side of (8.6) and in order to obtain the 
negative root, to substitute the value of © corresponding to the points 
Oo and 05 of the trajectory (see Fig. 5). Also taking into account the 
fact that tan @ < 0 at the points 2 and Oo of the trajectory and 
tan 9 > 0 at the points O5 and Oo, we can rewrite (8.5) for tan 0 
thus: 


r & 
(8.7) -tan@ = Ty (xtan 37 tan 1) 


On the basis of (7.16), we can write the equality 
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o = dy @ 
(8.8) X, tan eo tan a, = tan 6, + hy,cot 3 


which is valid for both ® and for ©'. Using it, let us rewrite (8.7) in the 
final form 


1 Basece 
(8.9) tan@ = - an a (in 0, + h,cot 4 


If it is necessary to determine the direction of the velocity at the 
point of intersection of an elliptical trajectory with a sphere of radius 
Ty then 


— ® 
(8.10) tan = - — (tan a, + h,cot 4 
1+ hy 2 


can be used in all cases. However, it is here necessary to take the value 
of ® determined from (7.35) for the points O, and Op» and the value of 


@' determined from (7.36) for the points Oo and Oo ; 


Evidently this derivation is valid for all three kinds of trajectories, 
elliptical, parabolic, and hyperbolic. 


9. FLIGHT TIME 


The determination of the flight time is related to the integration of 
the differential dependence expressing the law of the motion of the point 
along its trajectory. 


1. For a radial trajectory, the dependence 


dr 


= ty 
alt 


can be taken as such a differential dependence: hence there follows 


(9.1) dt = &t:— 
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Expressing v according to (5.2) interms of r, we obtain 


(9.2) ves -a 


b 
r 
where 


Zig 27g 
(9.3) b= 2T 95 a= yr, = ky) = 


, @ - &) 


In order to simplify the operation of integrating the dependence (9.1) 
it is convenient to transform from the variable r to the variable v 
rather than from v to r by means of the relations: 


et ee 


(9.4) Fe 
(a + v*)* 

which are easily obtained from (9.2). As a result, (9.1) can be transformed 
to the following: 


(a + v2) 


Let O, denote the point of maximum distance from the earth's 
surface on a radial trajectory by analogy with the corresponding point of 
an elliptical trajectory. 


For radial trajectories for which the point Op lies at a finite range, 
the flight time from the point O, to the point Oo can be measured by 
two methods: as the flight time from the point O, tothe point Oo along 
the shortest path, and as the flight time from the point O, to the point 
OF and from the point Of to the point Oo. Let Tj) denote the first 
flight time and T,9 the second. Then 


0 V 
2b dv 2 2b dv 
(9.6) ee Rae a7 eee 
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and 


V2 
ao ceetey for Vi 7? Vo 
12 9 2 
V1 (a + v*%) 
(9.7) 
V2 
Tio = aoe Nes for vy < Vv, 


2 
Vz (a + v’) 


since the velocity diminishes from v, to zero for motion from the point 
O, to the point Of and increases from zero to vo for motion from the 
point Or to Oo. 

For radial trajectories for which the point Of, is at infinity, the 
flight time between two points can be measured by only one method, as 
the flight time from the point O, to the point Oo along the shortest 
path; therefore, it will be determined by (9.7). 


The integral 


dv 
(9.8) I = 2b f—*<. 
(a + v2)" 


figures in (9.6) and (9.7) and it can be represented in the form: 


I {] 
es to 69 ho 
uo om 
—X, 
ae 
o 6 6hF IK 
= < 
bho es) 
+ 
7 tO |e 
“ 
a 
ae 
3 a a 
] + _—_ 
ee a“ 
Is 
—~, Si et 
LJ 
ao) 
a 
ww 
ee 


or 


(9.9) I 


I 
- 1c 
es 
— 
+ Je 
< 
te 
+ 
+ [= 
a 
a“ 
te 
| eee 
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For radial trajectories in which the point Oy lies at a finite range, 
k < 1; therefore, in conformity with (9.3), a > 0 always. Hence, the 
integral (9.9) can be reduced for them to the form 


2 
atv a a 


b 
(9.10) I= — [3 + —— avetan v2 
Moreover, in conformity with (6.2) and (9.3), the quantity a for these 
trajectories can be represented thus: 


on r 
0 1 r 
(9.11) a= ; or = ———— = 


Using (9.11), (9.3) and (9.4), let us write (9.10) as follows: 


r 
max k k 
: = EE + : 
(9.12) I = E J arctan ,/ | 


Now (9.6) and (9.7) can be reduced to the form: 


r k k 
max 1 1 
(9.13) Ty9 - k, + a ie a7 arctan Se oo. mare 
1 1 1 
r k Ko 
+ a ky + a/ = arctan ~/ 7 k 
Vo aS tS 
r k k 
' = max 1 1 
(9.14) Tio = + (., +a/ 1_k arctan ~/ 2 Ky 
r k k 
max y Meas. ey somes 
= (, + i kK arctan i - =| 


1 
where the plus sign in the last expression corresponds to the case 
Vy 7” Vg and the minus sign to the case Vi < Vo; 
For radial trajectories for which the point Of is at infinity, k =1 


and, therefore, a = 0 always, in conformity with (9.3). In the case 
‘a = 0, the integral (9.8) equals 
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2 
(9.15) I= - 4 


In the case a < 0, the expression (9.9) can be reduced to the form: 


iS ee (eee ne 
eive + a 2~./-a ve+a 
a | {se - Vi) p awe -_ 
ah + a 2./-a Vv - af v + a/-2 
or 
(9.16) I= 2 = + ~~ pm YY? 
ve+a 2./-a v + a/f-a 


Using relations (9.3) and (9.4) as in the previous case, let us write 


(9.16) as 
1 +2, {© in Vk-1 c - 


re 


r 
(9.17) I= = 


after which the flight time between two points of radial trajectories for 
which k = 1 is expressed as 


' v2 


and the flight time between two points of radial trajectories for which 
k > 1 is determined by the following expression: 
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r k k i 
y) 1 1 2 2 
9.19) T.. = +/—< i 
( 12 vo k.-1 27% k.- i ™ 
2 y) kK, 
: 7 eo 
‘Qo 7 
ky 3 
ee! : k, ++ “1 In = al 
+ 
k, - l ; 


The plus sign in the written expressions corresponds to the case 
Vv, 7 Vg and the minus sign to the case Vy < Vo- 


2. For the circular trajectory it is convenient to take the dependence 
(3.3) as the differential dependence from which the flight time can be 
determined simply. Rewriting it in the form 


NO 


Er 


(9.20) dt = = dy 


we obtain a very simple expression for the flight time between two points 
of a circular trajectory 


5s ee ee, ee 
(9-21) Ty, = GM. - 1) =F S42 = 7 Me 


(9.22) T, == =>. 


3. For an elliptical trajectory it is also convenient to take the 
dependence (9.20) as the differential dependence from which the flight 
time can be simply determined. However, it seems to be more complicated 
to integrate than in the previous case, since the dependence relating r 
and 1) is more complex for the elliptical trajectory. This dependence is 
expressed most simply by means of (3.26). Using it, we can express r 
in (9.20) in terms of 4 and we can integrate with respect to the variable 
nN, OY we can express dy interms of r and dr and integrate with 
respect to r. 
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Let us first proceed according to the second path. Solving (3.26) for 
7) we obtain 


1 p 
(9.23) 7 = Ny = arccos = ( - 3 


from which 


(9.24) dn = +—— ar 


where 


(9.25) Aw) = .A(L - e*)r* + Qpr - p? 


For elliptical trajectories (1 - e*) > 0; hence, it is convenient to 
transform (9.25) as follows: 


Jf A(r) 


Substituting (9.26) into (9.20), we obtain 


0” 
V1 - ec rodr 


= 4 
(9.28) dt = + = 


It is evidently most convenient to integrate the obtained dependence 
with respect to the variable x. To do this, let us express r and dr in 
terms of x. From (9.27) we obtain 


(9.29) r= eae + ex) 


(9.30) dr = aC, 


3] 
- eo 


Using the last relations, we can rewrite (9.28) thus: 
2 
p 1 + ex 


ao(1 - @2)8/? 1 - x? 


(9.31) dt = + dx 


Integration of (9.31) is not difficult. It is necessary only to establish 
the limits of integration first. We can rewrite (9.26) as 


r % 
(9.32) an + 5 

It follows from (9.32) that the sign in front of the square root when 
integrating (9.31) must be chosen with consideration of the sign of the 
derivative dr/dn on that section of the trajectory over which the inte- 
gration is made. Let us call the part of the trajectory on which dr/dyn > 0 
the ascending branch and the other part of the trajectory on which 
dr/dyn < 0 the descending branch. Evidently the points at which 
dr/dn = 0 are the vertices of the elliptical trajectory O, and Oy 
and they separate the ascending and descending branches of the trajectory. 
It follows from (9.32) that the values of x 


(9.33) %* = +1 and x = -l 


correspond to the above-mentioned points. Substituting them into (9.29), 
we obtain the values of r corresponding to the vertices of the trajectory: 


(9.34) r= and r = 


which agree with the values of (6.4) described earlier. Here x = +1 
corresponds to r = Tryax i.e., the point Of, and x = ~-1 corresponds 
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to r = Tin, i.e., the point Op. Evidently the plus sign must be taken 
on the ascending branch of the trajectory between the points Op and OF 
and the minus sign on the descending branch between the same points. 

As in the case of the radial trajectories, we denote the time measured 
between the two flight points O, and O» which are on different sides of 
Of by T, 9 and define it as the flight time from the point O, to the 
point Op andfrom Of to O9. Then from (9.31) follows: 


x 
(9.35) Tyo = p* i (1_ + ex)dx r (1_+ ex)dx 
’ 12 — ~ 
o(L - e4) x, vi - x7 1 V1 - x* 
where 
(Ql - er, - p (l - er, - p 
(9.36) Xx, = pe? Xo = Se 


As before, let Ti denote the time measured along the shortest 
path between the points O,; and Os» on one side of the point Of. On the 
basis of (9.31), we can then write 


2 X9 
+ 
(9.87) Ty = —_—-__ (i + ex) ay for ry <r 
o(1 - 0283/2 xX fl - x2 : 
1 
- p* 72 1 + ex 
(9.37') Ty, = Hf M+ 0) ay for rj > Yo 


o(1 - 9238/2 X1 J/1 - x2 


where X, and x, are also determined by (9.36); the expression (9.37) 
corresponds to the case when both points are located on the ascending 
branch of the trajectory and (9.37') to the case when they are on the 
descending branch. 


The time for a complete revolution on an elliptical trajectory is 


‘) 


—_ 


+1 -1 
(9.38) 1 = p TE (1_+ ex)dx _ (1_ + ex)dx 
o(1 - e*) -1 Vi - x? +1 V1 - x? 
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Using the fact that 


(9.39) f 28 ax = f—K~-, (ges 


JL - x? J1 - x2 afi - x? 


arcsin x - e 1 - x? 
we can represent the flight time for a complete revolution as follows: 


2 


p 
O40). 1.223 
oa - e2)3/2 


after which the expressions for Tyo and Tio can be represented as: 


(9.41) Tyo = Ty, ( - aresin x; + evi - x4) 
+ t - arcsin Xo + et - 3) 
(9.42) Tio = + Ty, (GE - arcsin x; + evi - x4] 
_ : - arcsin Xo + e/i = x8)| 


where 


T 2 
Cc Pp 
(9.43) T,. 2. —— 2 SS 

27 o(L - 92 3/2 


is the average time corresponding to the rotation of a point through one 
radian around the center of attraction.* 


*Briefly, the mean radian time; in conformity with this we append the subscript 
r thereto. 
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Using expressions for 9 and a in the form: 


(9.44) 0 = /Top; a = 


Before reducing (9.41) and (9.42) to final form, let us integrate the 
dependence (9.20) with respect to the variable 7 by expressing r in 
terms of 7 by using (3.26). As a result we obtain: 


ys 
p 
(9.46) dt = — a : 
k - e cos(7 - "p)| 
or 
2 
(9.47) dt = = et = 
E + e cos(7 - 7 } 
F 
where 


(9.48) Tp =n, +7 


is an angle which is more useful than the angle Np in astronomy. 


Let a denote an angle measured from the vertex of the trajectory 
Op in the motion direction; then we can write 


(9.49) a 


ll 
3 
' 
3 


t 
F 
and 


TO ae eee 
(1 + e cos a)? 


Dividing numerator and denominator on the right side by cos* a/ 2. 
which is valid for all values of a@ with the exception of 


(9.51) a@ = Nn 


where N is any integer, we obtain 


s 
1 a($) 


22 2% ( 2 .) ( e) 
1 cos 9 COS” 5 _ 2p? 1 + tan” 5 ditan > 
~  ¢ | a 2 
-— te 1 + tan® = + ef - tan” 4 
cos? = cos? — a ? 
2 2 
*s)een§) gt _(+ mt) aon) 
: 2p* (1 + tan )a(tan 9 ; an 1 + tan d{tan 9 
2 2/2 oo + e)2 l-e, 2/2 
Ja+esa eytan® 5] ( ) 1+ >> tan’ 5 
Multiplying and dividing the right side by | ; : = > 0, let us 
rewrite the expression for dt: 
l-e ,1l1-e 2 2) af 1-e 2) 
ie 2p* & fe Tee BW iee. 8 
- 3/2 l1-e, 9a/2 
o( + e (+=) fn + Pe tan 
lt+e 
Then, using the notation 
7 - nN 
€ 1-e a 1-e 
—_— = — t 5 
pe “ERS ite * 2 Ite 2 


we 


1l1-e 


(Ces 


+ (1 + tan” <)| a(tan 


can represent the expression for dt as follows: 


€ 


ee 2p* /fi+el\l+e 2 
7 2V1-e 2«€)\4 
o(jl1 - e) 1 + tan 9 
ze E € 
— t is 
2p° l+e cre atm) alten 5) 
= ~ 2 
- € € 
gee P= (1 + tan? 5) 1 + tan? — 
2p" lte 2e 9 € 1 1 
x 1 "Tee 8 boo. Sate 
v(1 - e*) . 
p? lt+e e 
= a — Tee Ut cos €)de + de 
o(1 - e”) 
p? 
= (dé - e cos € dé) 
3 
a(1 - e%) /2 
Finally, using (9.43), we obtain the final expression for dt: 
(9.53) dt = T, (d€ - e cos € dé) 


Integration of the dependence (9.53) is not difficult; it is only 
necessary to establish the limits of integration. In contrast to the 
preceding case, the integrand has no singularities. Hence, to establish 
the limits of integration it is sufficient to clarify which values of «€ 


correspond to the points O 
the value of € 1 determin 


€ 
Li... fee l 
(9.54) tan 2 ar tan : 


and O,. Returning to (9.52), we obtain that 
by the expression 


corresponds to the point O,, and the value of € 9 determined by the 


expression 


D0 


(9.55) 


corresponds to the point Oo, where, in conformity with (9.48) and (4.15), 
the magnitude of us is determined by the formula 


(Uor - 1)sin 7 + tan @ cos 7 
(56). Stang =. 
az28) Ip (Upr - l1)cos 7 - tan 0 sin 7 


Integrating (9.53) from the point O; to the point Oo; we obtain an 
expression for the tlight time between the points in the following form: 


(9.57) Tyo = a & - e sin €9) - (€, - e sin €y)| 


where €1 and €, are defined by (9.54) and (9.55); but (9.57) remains 
the same for any location of the points O, and Oo. 


Considering the parameters r., n,, v,., 9, prescribed at the point 
O, and the parameter ro atthe point Oo, we reduce the formulas 
obtained for the flight time to a form convenient for computations. 


Substituting (4.12) and (4.14) into (9.36), we obtain: 


- 2k 3 
1 1 2k, 


(9.58) xX, = eee. z 


after which (9.41) and (9.42) can be represented as follows: 


(9.59) Ty. = 7,16 - By + ¥;) + € - By + %) 


(9.60) Ty. = «1,|(3 - By + 11) 7 (3 a v2) 


The plus sign in the last formula corresponds to the case ry < %, and 
the minus to the case Yr; 7 Yo and in both formulas 


a3 


1 - 2k 


i { 1 a < Tv 

B, = arcsin x, = arcsin —— (- ie B, = ; 
1 - 2k : 

B. = arcsin xX». = arcesin orn : ; 


(9.61) vy, =e /1 xy 4k, (1 k)sin a, e cos B, 
Yo =e, /1 Xo / 4k, (1 k,)sin 9 


ae ees ee ee 
— 201 - ky) 2(L - ky) | 


| 
© 
OQ 
fo) 
n 
D 
iw) 


When using (9.57) it is convenient to transform from the angles 7 
to the angles f, defined by (7.3), in (9.54) and (9.55). Since (9.54) and 
(9.55) are valid for any case of measuring the angles 7, then they are 
valid also for the angle f which is the angle 7 measured from the 
radius vector to the point O,. Replacing the angle 7 by the angle f, 
we obtain 


@' 
eee ae ae “4 
2 1+e 2 
(9.62) b- 6 
a ee pe ae ee 
2 1+e 2 


where in conformity with the expression (7.32) 


g 
2k tan 1 


' = 
(9.63) -tan®, = ; 


1 


4. For the hyperbolic trajectory the expression for the flight time 
can be written by using the method used for the elliptical trajectory. 
Hlowever, in transforming the expression for /A(r) it is necessary to 
take into account that e > 1 and to transform it thus: 


/A(r) 


= 1)r? + 2pr - p@ 


(e 


- l 
(9.64) dt = + a 


where, just as in the previous case, we can put 


2 
- + 
(9.65) 2 = S— FIP. 
pe 


It follows from (9.65) that 


(9.66) r= se - 1) 
e“ - 1 
and 
e 
(9.67) dr = —Pe dz ; 
e* - 1 


(9.68) dt = + ——— — > a 
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The value z = 1 corresponds to the vertex of the trajectory Oy 


p : 
for which r = ee hence, reasoning as in the preceding case, we 


obtain an expression for the flight time measured between two points of a 
hyperbolic trajectory on different sides of the vertex OF in the form: 


1 
(9.69) ty = — a | —_ Shel 


o(e# - 1) ay; Ce 


and for the flight time measured between points on one side of the vertex 
O, of the trajectory, in the following form: 


p* -2 ez - i 
(9.70) Tyg = + ————7ZR f —<={ «wz 
oe? - 1%" zy vo% - 1 


The plus sign corresponds to the case r, < r, and the minus sign to 
the case r, >. ro (it is assumed that the motion along the trajectory is 
from the point O; to the point Oo). 


Using the fact that 


[ Fain e fiat. JT=- er/z* - 1 - arccosh z 


eJ/z* - 1 - Inz + V2" - 1) 
we represent the expression for Ty and D5 in the final form: 


Ty0° = T,. ( 7 - 1 - arecosh 2) + (« Ver - 1 - arccosh 2] 


(9.71) 


toto 


ae = +£ T,](es Z5 - l - arecosh 9) - (e73 - ] « arecosh 4) 


or 


(9.72) 
' = f2 2 
-<e@4/z, - 1 - In(z, + | 
1 1 1 
where 
T. = pe 
ro 3 
ol1 - e4| [2 
2 
(e - lyr, + Pp 2k, - I 
(9.74) Z4 = a = a 


ae a 
(e Ir, + p 2k, 1 
Zo SSS SS 
pe e 


Just as before, going along the second path and integrating the 
dependence (9.46) with respect to 7, taking into account the fact that 
e > 1 for a hyperbolic trajectory, we obtain instead of (9.53) 


(9.75) dt = T. cot c : id 


cos” € 


where € is defined by the expression 


€ G 
Pf). aD. e+i1 2 /e+l 2 
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which is analogous to (9.52). Integrating (9.75) from the point O01 to Oo, 
let us write the expression for the flight time between these points as 
follows: 


(9.77) Ty = T,<je tan €5 - In an\ 7 9 


€ 
T 1 
- fe tan €4 - In tan(Z + 2} 


where ey and €, are defined by the expressions 


c t € ® = @! 
(9.78) tan oa ( tan =), tan — = tan ——— 


which are analogous to (9.62), and the quantity On by (9.63). 


5. For a parabolic trajectory the expression for the flight time can 
be written also by following the method we used for the elliptical and 
hyperbolic trajectories. It is necessary to take account only of those 
simplifications, which from the very beginning, result from the equality 
e = 1 for this case. 


Thus, on the basis of (9.20), (9.24), (9.25), we can write for the 
e = 1 case: 


(9.79) dt = + 


J2pr - p* 


Exactly as in the preceding cases, the radicand vanishes at the vertex 
of the trajectory O, for which 


p p 
(9.80) Yr = Ymin = a 


1 +e 


Hence, by analogy we can write the expression for the flight time 
measured between two points of a parabolic trajectory on either side of 


the vertex OF as: 


Pmin 7) 


d 
- f = > if — 
ry a/2pr - p min -V 2PY - p 


oe ke) 


(9.81) T 


12 


and for the flight time measured between points on the same side of the 
vertex OF of the trajectory as: 


(9.82) Tye = 


where the plus sign corresponds to the cuse ri < To and the minus sign 
to the case ry, > To. 


Using the fact that 


y r_dr oot, eS ee ae 


o if = o 
/opr - p2 *P /2pr - p* 


1 [2 9.3/2 : 91/2 
apo F (2pr - p*) 2p“(2pr - p“) 


2 


1 5 _P (2 ) 2r 
gg F + P)fepr - P= 35\h tt [> 


let us represent the expressions for T,. and Tj» thus: 


2 [2x r or 
p ( 2 2 1 1 
(9.83) 12 aa | 7 + , 5 1 + (- + 7 " | 


2/r or r or 
; p 2 2 ( 1 1 
= —— —= 1 as — + 1 —_ 1 
(9.84) Ty) = + ale +1) : : : 


where 


HK 
| 


(9.85) p = 2r cos?6 ; 0 = rv cos 8 


Determining the expression for the flight time by the second means, 
i.e., by integrating (9.50) for the case e = 1, we obtain in place of (9.57) 


Q! Xo 
T _ p* ia da eS p* 7 
i200 ¢ (a +ecosa)? 9% & ey 
1 1 2 
re. 
OQ! d - 
Be il 2 
~ 20 a 9 & 
a, cos 9 COs” 5 
2 2 
_P ( 2 2) ( 7 
ao f 1 + tan® 5} d\tan > 
“a 
or 
(9.86) Tyo 9G (tan > + 3 tan =) 7 (tan Pa tan + 
where 
(9.87) a, = - On; d, = &- OF 


in which ®, is defined by (9.63) and 


® = No - 14° 


10. IDENTICAL VELOCITY TRAJECTORIES. 
MAXIMUM RANGE TRAJECTORY 


Let there be prescribed the coordinates of the point O1 and the 
magnitude of the velocity Vv, at this point, i.e., the quantities 


(10.1) Yy> Ny> V1 


60 


Then, by prescribing various launch angles 04 we obtain a family of 
trajectories, to each of which corresponds its value of the angle o,- 
This family of trajectories, emerging from the same point O, and 
corresponding to the same value of the velocity at this point, will, for 
brevity, be called the family of identical-velocity trajectories. Certain 
properties of identical-velocity trajectories will be examined in this 
section. 


If we are given a certain sphere of radius r, measured from the 
center of the earth, then the question can be posed of the determination of 
points of intersection between the trajectories of the family we are 
interested in and the above-mentioned sphere. 


Let us analyze the following cases separately in solving this problem. 


1. We will first examine the case when the radius ro of the given 
sphere is less than the radius r, and all the trajectories of the family 
intersect the sphere. Evidently this case corresponds to absolute 
incidence of the moving point of interest onto the prescribed sphere when 
a certain constant velocity is imparted (to the point) at the point O, 
located on the sphere. Since the family of parabolic and hyperbolic 
trajectories passing through the point O; will always contain a tra- 
jectory going from Oj, to infinity, then only elliptic trajectories can be 
the trajectories which will absolutely intersect the prescribed sphere. 
Let us examine which conditions the family of elliptical identical-velocity 
trajectories must satisfy in order that all its trajectories should intersect 
the assigned sphere. 


If all the elliptical trajectories of the family intersect the prescribed 
sphere, then the angular range ® of the point of intersection of the 
trajectory with the sphere can be put into correspondence with each of 
them. Hence, of the two points of intersection of the ellipse with the 
sphere which exist in the general case and which were denoted for the 
case under consideration by Oo and Oo on Fig. 5, let us put the angular 
range of the point Oo, defined by (7.35), into correspondence with each 
trajectory of the family. According to the prescribed parameters (10.1), 
which define the family, and the prescribed angle 01, which determines 
the chosen trajectory of the family, this formula always permits finding 
the angular range & corresponding to the chosen trajectory. This same 
formula permits a sufficiently simple determination of the condition 
which the family of identical-velocity elliptical trajectories must 
satisfy in order for all its trajectories to intersect the given sphere. 
Formula (7.35) shows that there will be no intersection of the trajectory 
with the given sphere when the radicand is negative. Vanishing of the 
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radicand means coincidence of two points of intersection at one point of 
tangency of the trajectory with the sphere. To clarify the conditions 
under which the radicand in (7.35) vanishes or becomes negative, let us 
substitute therein (7.17) for X,. Asa result, the desired condition of the 
presence of a point of intersection can be represented as follows: 


2 eS eee, ee " 
(10.2) soothe le + Ay - 2-h,} > 0 


and the condition of absence of a point of intersection as: 


h 
2 i - (1 : 
(10.3) wae ea < 0 


In the case under consideration, Ty > Lo 


1 + — > 0; consequently, the inequality (10.3) can be reduced to 


and, therefore, h, > 0 and 


(10.4) tan*6, < 


from which it follows that under the conditions 


(10.5) 2+ hy < ork < *— 
1 2+ hy 
the inequality (10.4) is impossible for any value of ¢, since the left side 


is always a positive quantity and the right side negative. But this in turn 
means that under this condition there are no trajectories among the tamily 
under consideration which would intersect the sphere. 


Hence, the first case we have examined of a family of identical- 
velocity trajectories must be characterized by the conditions: 


(10.6) h, > 0; k, < 


1 
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Therefore, a definite flight range, computed by means of formula (7.35), 
corresponds to each of the trajectories of this family. Evidently zero 
angular ranges will correspond to the two trajectories of the family which 
correspond to the angles 6, = 1/2 and 6, = - 1/2 and, therefore, 
degenerate into radial trajectories. It is also clear that upon the transi- 
tion from a radial trajectory corresponding to the angle 9, = 1/2 to 
elliptical trajectories for which the angle ¥, gradually diminishes while 
remaining close to 7,/2, we will obtain an increase in the angular ranges 
correspondiny to these trajectories. In exactly the same manner, when 
making the transition from a radial trajectory corresponding to the angle 
0; = - 7/2 to elliptical trajectories for which the angle 4, gradually 
increases while remaining close to - 7/2, we will also obtain an increase 
in the angular ranges corresponding to these trajectories. Consequently, 
it is natural to expect that an angle oy is found between the angles 

- 7/2 and 7/2 to which the maximum angular range will correspond. 


The problem of determining this angle mathematically requires 
determining the value of the angle 9, for prescribed r,, vy and ro 
for which the function 


(10.7) & = O(vz, 9, Ty lo) 


takes on the maximum value. The solution of this mathematical problem 
reduces, as is known, to looking for the roots of the equation 


o®(6,) 


(10.8) 50, = 0 


To obtain the necessary partial derivative of the angular range it is 
simplest to return to (7.16) and, assuming k, and h, to be constants 
and ® to be a function of only 9, to differentiate it with respect to @,. 


As a result we obtain 


d ce) 
ae — @ ; 
a > o : 2tan — ad 2 tan | 2 tan l ge - 
—— tan” — a ee a ee ee 
1 ® ® J 
oe 2 2c0s* os oe cos*6, 2cos*— 1 


where on the basis of (7.17) 
ox 2 sin 0, 

2 
1 k,cos a, 
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Solving the obtained equation with respect to 0®/ 06, , we write 


-1 
= tan 9 tan = 
I® : 2 tan o4 tant 2 2 tan 9 1 7 9 
00, 2 2 2 OE ee 
1 k,cos 04 cos a, cos 9 cos 5 
or 
® co ® 
e eee panes @ pare 
36 cian eae 1 tan > a l 
00, k ® 


Using (8.8) and (8.10), let us rewrite the obtained expression as 
follows: 


® 
—_— 4 rc) 
‘ea ab (tan a, tan , K,) sin 
80, r £) 2 
1 (tan 0, + hycot ; k,cos 0, 
or 
tan 0, ¢t = k in ® 
06 an 9, tan aoe :) sin 
(10.10) 3e. 7 7 ; 
1 (1 + h,)tan a4 k,cos a, 


Hence, it follows that the derivative 90/ 06, vanishes for ® # 0 only 
when 


Qn 
(10.11) Kim = tan Oam tan cS 
This will evidently also be the desired relation between the trajectory 
parameters which will guarantee maximum range. Later, just as has 
been done above, we shall ascribe the lower subscript to the right, m, 
to the trajectory parameters which guarantee maximum range, and the 
trajectory itself shall be called the trajectory of maximum or greatest 
range. 


G+ 


The desired thrust angle can be determined from (10.11) and we shall 
call it the angle of thrust of the maximum or greatest range; it will be 


denoted, as above, by Om! 


Pn 
‘ 0 = nae 
(10.12) tan ies Ky moot , 


However, it is still impossible to determine the angle Om by means of 
the given parameters (10.1) by using (10.12); to do this it is necessary to 


find ¢,,. 


In order to find %,,, let us substitute (10.12) into (7.16). We obtain 


b om 

2 ee _m _—m _Ff - ‘ 

X imt22 2kim cot - tan - hy 0 
where 
2 
, 1 _ 1 + tan om _ 
Im * ; - (2 + hy) = k - (2 + hy) 
} 
1 = 9 *m 
-(- -2- i) +% cot” — . 
oe i - 


Hence, the equation to determine %,, is reduced to the form 


2 2m h 


2 


(10.13) te Oe h, tan 


Kim 


= Kim 


from which 


@ 
(10.14) tan = 


65 


Substituting (10.14) into (10.12) we obtain the formula to determine 1, 
for the given parameters (10.1) in the form 


(10.15) tan, = 


Formula (10.14) permits a particularly graphic answer to the 


question: Can a point with a velocity v, at an altitude h, move along 
an elliptical trajectory enclosing a sphere with the radius ro described 
from the center of the earth? 


Evidently, an elliptical trajectory will enclose this sphere when 


1 


(10.16) ky, > ———, 
mm 2+ hy 


i.e., when the right side of (10.14) is imaginary and, therefore, the point 
of intersection of the trajectory with the sphere is lacking. 

We shall turn below to the consideration of a case where the maximum- 
range trajectory degenerates into a trajectory enclosing the prescribed 


sphere. 


2. Thus, let us examine the case when the family of identical-velocity 
trajectories is characterized by the following conditions: 


(10.17) h, > 0; k, > onli 
2+ h, 


As follows from (10.4), points of intersection with the prescribed sphere 
are lacking in this case for trajectories to which correspond values of 
04 satisfying the inequality 


(10.18) J] < 4% 


where 


(10.19) 93, = arctan 
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«And points of intersection with the prescribed sphere are natural for 
trajectories to which correspond values of o4 satisfying the inequality 


(10.20) 1°, > O14, 


Each of the latter trajectories can be put in correspondence with the 
angular range & of the point of intersection which is foremost among the 
forward-lying points of intersection. This angular range, as has been 
shown in section 7, is determined by (7.35) for all kinds of trajectories in 


the case hy > 0. 


For all kinds of trajectories a zero angular range will correspond to 
the angle 6, = - 7/2. As the angle 0, increases from the value - 7/2, 
an increase in the angular range can be expected as long as the angle 04 
is not equal to - ,,, for which the point of intersection transforms into 
a point of tangency.* If a1 = - 04, is substituted into (7.35), we obtain 
the following formula for the angular range of the point of tangency: 


® 
(10.21) tan =. eee 
aes 


also valid for all kinds of trajectories. As the angle 9, increases 
further from - @,, to @ ,? we shall obtain trajectories enclosing the 
prescribed sphere and not intersecting it. As the angle o4 increases 
from 1 to 7,2 for elliptical trajectories, we will obtain the forward- 
lying points of intersection, where the angular range to the foremost of 
these will be determined, as before, by (7.35). At the point of tangency 
which is characterized by the value of the angle i = 94h this angular 


range will be determined by the expression 


@ 
1 
(10.22) tan _2k SS ta eS 
2 tan 4k 


*The subscript 'k" has been ascribed to this angle to emphasize that it 
corresponds to the tangency case (Russian word ''kasaniya'' for tangency). 
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As the angle 6, increases from the value Qa? a decrease in the 
angular range to zero for the angle eT = 1/2 can be expected. As the 
angle o, increases from 91, to m/2 for the parabolic and hyperbolic 
trajectories, in conformity with section 7, we will obtain only the back- 
ward-lying points of intersection of the trajectory with the prescribed 
sphere which will not be reached by the point under consideration moving 
from the point O,. Hence, it has no meaning to set these trajectories in 
correspondence with some angular range. It is more natural to consider 
these trajectories as not intersecting the prescribed sphere since they 
recede from it to infinity. 


Hence, in the case under consideration, as the angle @, varies 
between - 7/2 and - 9,, the angular range varies between zero and 
the angle ,),, which is determined by (10.21) and satisfies the condition 

b = Tf 
[eux] = * 


As the angle 6, changesfrom - 9,, to 4},, the trajectories of the 
family do not intersect the prescribed sphere. As the angle 6, changes 
from 64, to t/2 the elliptical trajectories intersect the prescribed 
sphere and the angular range changes from Po. to zero. Here 5, is 
determined by (10.22) and satisfies the condition ®5, = 27 - %4,. As 
the angle 91 changes between 6, and 7/2 the parabolic and hyperbolic 
trajectories do not intersect the prescribed sphere in forward-lying 
points. 


The monotonicity of the change in the angular range in the mentioned 
intervals of variation of the angle 9, is easily discerned physically and 
can be established mathematically by an appropriate analysis of formula 
(10.10) for the partial derivative 9%/00, . 


3. Third, let us examine the case when the radius re of the 
prescribed sphere is larger than the radius ry and all the trajectories 
of the family do not intersect the sphere. Evidently only elliptical 
trajectories can be the trajectories that do not intersect the sphere in 
this case. In order to determine what conditions must be satisfied by the 
family of identical-velocity clliptical trajectories which do not reach the 
prescribed sphere, let us return to inequality (10.2) which is the condition 
for the presence of points of intersection of the trajectories of the tamily 
with the prescribed sphere. Let us rewrite this inequality thus: 


h, ee 
(10.23) tan®@,\1 + 7} > (-hyl - 2 - hy 
‘1 4 
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and taking into account that -1 < hy < 0 always for r, < Yo, letus 
examine two cases separately. We first assume that 


h 
(10.24) 1 + at > 0, i.e., = rer or - hy < k 
k k h 1 
1 1 1 


Then the inequality (10.23) can be rewritten as 


P (hy) /y : 
10.25) tan“@, > ————|(— - J 
( ) 1 b, ik 2 ) 
1 + 2) 
1 

Then, we assume that 

h 
(10.26) 1 + 2 < 0, i.e., Ps a or -hy > ky 

Ky Koh 

1 


Inequality (10.23) can then be rewritten as follows: 


Y /1 . 
(10:27): tane0,, <6 | 2 Ty 


It follows from the written inequalities that the condition of absence 
of intersections of trajectories of the family with the given sphere must 
be the inequalities 


(10.28) oe ieee ae oe hy 


for compliance with which the left side of inequality (10.27) must always 
be a positive quantity and the right side, negative, which indeed means 
the absence of trajectories intersecting the given sphere. 
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Let us determine whether the second inequality of (10.28) is a 
consequence of the first or not. To do this, let us verify the validity of 
the inequality 


For h, < 0 it easily reduces to the always valid inequality 
- 2 


Hence, the inequalities (10.28) can be reduced to the one: 


k < -h, 


and the considered third case of the family of identical-velocity 
trajectories must be characterized by the following conditions: 


(10.29) -1 <h, <0; k, < -hy 


4. The last, fourth case of a family of identical-velocity trajectories 
which is characterized by the absolute presence of trajectories which 
intersect the given sphere under the condition h, < 0 will evidently 
hold upon compliance with the inequalities: 


1 1 e 
(10.30) erp k, > -hy, 
1 hy 


Here, if the inequalities 


1 m 
(10.31) me > ky > -h, 


2+ hy 


are satisfied, then the elliptical trajectories corresponding to these 
inequalities and angles eT included within the interval 


(10.32) 4 < O44 
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where 


(10.33) OK = arctan 


do not intersect the given sphere, and elliptical trajectories corresponding 
to the same inequalities and angles 4 included in the interval 


(10.34) |0,| > 64, 


intersect the given sphere at two forward-lying points. The angular range 
to the first of the points is hence determined by formula (7.36) and the 
angular range to the second, by (7.35). Evidently, as in the case con- 
sidered earlier, the angular range to the first point of intersection will 
change from zero to #1, asthe angle 6, changesfrom - 1/2 to 

- 94, and will change from 5, to zero as the angle 9, changes from 
Qk to 7/2. Because of the agreement of the expressions for O4;, in 

the present and previous cases, the expressions for Pik and Po. will 
also agree with the expressions (10.21) and (10.22) obtained in the 
previous case. 


When the inequality 


(10.35) k, > — 
2+ h, 


is satisfied, all the trajectories intersect the given sphere, where the 
elliptical trajectories intersect the sphere in two forward-lying points, 
the angular ranges to which are determined by the same formulas as 
above and the hyperbolic and parabolic trajectories intersect the sphere 
in one forward-lying point, the angular range to which is determined by 
(7.36). 


The analysis of the possible cases of the intersection of a family of 
identical-velocity trajectories with a given sphere can thereby be 
terminated. In conclusion, however, it is interesting to examine the 
question of the envelope of the family of trajectories under consideration. 
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Let us recall* that a line tangent to different curves of a family at 
all points is called the envelope of the family of curves. If the family is 
prescribed by the equation 


(10.36) F(x, y, a) = 0 


where x and y are coordinates of points of the curves and the quantity 
a is a parameter of the family, then the envelope of the family of curves 
is determined by the equations 

uF 


(10.37) F(x, y, a) = 0; a F, (x, y, a) = 0 


In our case, it is convenient to take the trajectory equation in the 
form (7.5) as the equation of the identical-velocity trajectories in which 
u and f are coordinates and the parameter must be the angle 0,. The 
rest of the parameters are constants, by definition of the family. 
Rewriting (7.5) in the form of (10.36), we obtain 


(10.38) Fu, f, o4) =u + U,cosf + Uosinf - Ug = 0 


where Up, Uy, Us are defined by (7.14). 


In order to write the equation of the envelope, let us find the 
derivative Fy = oF /006,: 


guy oU wy 
og 


(10.39) F,(u, f, 94) = 
1 


Substituting the following expressions which are easily obtained from 
(7.14) 


(10.40) =— = —- = ———_.. = 
1 


“See, for example: V.T Smirnov: Kurs vysshei matematihi (Course in Higher 
Mathematics), vol. 2, Gostekhizdat, 1951, p. 40. 


2 


into (10.39), we obtain the desired equation in the following form: 


tan aT 1 tan o4 
(10.41) F,(u, f, 91) = ———z_ cos f + ———— sin f - 
kyrycos 0, r,cos o4 kyrycos*U, 


Equations (10.38) and (10.41) are indeed equations of the envelope. 
They can easily be reduced to a single equation if the parameter 6, is 
expressed from (10.41) as: 


k,sin f 


f 
: @ = le ee = ve 
(10.42) tan 6, io cost ky cot 9 


and then it is substituted into (10.38). To simplify the computations, it is 
convenient to represent (10.38) in the form: 


i wy tone & 1\ ef _ 
(10.43) (u - x) 2U sin > cos > - 2 (vu, = 1) sin = 0 


Substituting (10.42) into (7.14), we will have 


1 + 2 asic 


1 2 
Ug = 2kiry 
(10.44) ; 
k,cot > 
Us ee 
1 


After this, substituting (10.44) into (10.43), we obtain the desired equation 
in the following form 


ai ys ee ee 2 £ 
1) : 2k, cos 5 sin 5 + k,cos 9 : 2 sin 9 ae 
ry ry kyry Ty 
or 
k. - l 
10.45 ee ae er ae au Nea dine a 0 
(10.45a) |u ry r, 9 ty 9 
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Now transforming from the half to the full angles, we can write 


1 - ki 1-k, 
(10.45b) u = —>—— ( - =—— cos f 
or 
Po 
(10.46) r= 7 cos f 
) 

where 

. P, = ; ea. a 

0 0 
1 - Ke 1 + ky 


It follows from the expression for p, that the envelope exists only 
for elliptical trajectories when ky < 1. This result could easily have 
been foreseen; since it is clear from physical considerations that the 
families of parabolic and hyperbolic trajectories cannot have envelopes. 


The obtained equation for the envelope of elliptical trajectories shows 
that the envelope is also an ellipse whose focal axis is directed along the 
line AO, connecting the center of the earth A with the point O, which 
we have selected as the initial point for the family under consideration. 
Hence, as should have been expected, the enveloping ellipse is located 
symmetrically with respect to the line AO,. The semi-axes of the 
enveloping ellipse equal: 


+ ° 
me ee, oe 
0 =" =e 
1 - ee 2 1 ky 2ey 
0 
(10.48) : 
— ~ ay Ve 
0 Poo ~ *1 4-4 kK, 


As is known, the distance between the foci of an ellipse equals 
2a&9: For the enveloping ellipse 


(10.49) 2ae) = x, 
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and since one focus of the enveloping ellipse is at the center of the earth 
A, then, on the basis of (10.49), the second focus must be at the initial 
point of the family which we selected, O,- 


Different families of elliptical trajectories and, therefore, different 
enveloping ellipses whose foci will, however, always coincide with the 
points A and O,; will correspond to different values of the initial 
velocity V1: Hence, by varying V, we obtain a family of confocal 
enveloping ellipses whose semi-axes grow without limit as the velocity 
increases and approaches the parabolic velocity. Several such enveloping 
ellipses are shown in Fig. 7. 


The enveloping ellipse has an interesting property: none of the 
elliptical trajectories of the 

family under consideration emerges 
beyond its limits. Consequently, 
with its aid it is possible to 
estimate graphically the properties 
of the elliptical trajectories 

issuing from the selected initial 
point O, and corresponding to a 
given value of the initial velocity 


Va 


For example, by determining 
the angular range of the point of 
intersection of the enveloping 
ellipse with a sphere of radius 
r, described from the center of 

Fig. 7 the earth, we obtain an expression 
for the maximum angular range 
achieved by the family under consideration on the given sphere. To 
determine the mentioned angular range, it is sufficient to substitute the 
value u = 1 rp into (10.45b) and to solve the equation thus obtained with 
respect to the quantity f, which in the case under consideration will also 
be the desired angular range ©. As a result we obtain 


1 + k, 2k 
(10.50) cos #,, = 7711-7 7G 
te i-k 
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This formula does not disagree with formula (10.14) obtained earlier. 
Actually, if we return to the equation of the enveloping ellipse in the form 
(10.45a) rather than in the form (10.45b), we obtain the formula 


2m 


(10.51) tan —— 


which agrees completely with (10.14) and is valid both for To S$ Ty and 
for To > Ly. Inthe rp > r, case the above-mentioned formula yields 
the maximum angular range for the second point of intersection of the 
elliptical trajectory with the given sphere. 


11. IDENTICAL RANGE TRAJECTORIES. 
MINIMUM VELOCITY TRAJECTORY 


Let a certain initial and a certain final point of the motion O, and 
Og be assigned. Let us determine the magnitude v, and the direction 
U, of the initial velocity which would guarantee getting from the point 
O, to point Oo. 


Let us give the position of the points O, and Oy in polar coordi- 
nates 


(11.1) rj, f, = 0 and r,, f, = @. 


Then the condition for the trajectories under consideration to pass 


through the points O, and O, can be written in the form of (7.16) in 
which the quantities 


(A122). iy ty. @ 


should be considered known now and the quantities Vyze Uy unknown, 

Since there is only one equation for the determination of two unknowns, 
there is an infinite number of solutions. For example, it is possible to 
assign values of Uy arbitrarily and to determine the appropriate value 
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of v, for each. Hence, in this case also we obtain a family of trajectories, 
each of which is determined uniquely by the initial thrust angle 6. and 
which passes through the assigned points O, and Oy. For brevity, 
henceforth we shall designate these trajectories as identical-range 
trajectories. 


Solving (7.16) with respect to x1 and then with respect to k,, we 
obtain: 


i) — co) 
(11.3) X, = 2 tan Oy cot > + hy cot? 2 


and 
. 2 2 = ~ ye 2 
(11.4) Ky = sin 20, cot 9 + cos a, (2 + h, + h,cot :) ; 
and since 
2 
— “ae 
1 279 
then 
2 De a 
(11.5) vt = eo ta De loca. + fycot” )| 
r, |sin 2¢, co 5 cos“@, | 2. 1 1° , 


The obtained dependence 


(11.6) vy = vy (64) 


permits a determination of the magnitude of the velocity for a given 
direction, which is necessary to guarantee passage of the trajectory 
through the prescribed points O, and Oo. In particular, this dependence 
permits the determination of the minimum velocity possible for passage 
of the trajectory through the prescribed points. To do this it is sufficient 
to take the derivative of the function (11.6) with respect to 4, and to 
equate it to zero. The root of the equation 


dv, (?,) 
dé, 


(11.7) 


defines the angle 9,,, corresponding to the desired minimum velocity. 
Substituting this angle into (11.5), we obtain the magnitude of the minimum 
velocity. 


However, it is simpler to look for the maximum of 1/k, rather than 
for the minimum v, since the found maximum of 1/k, will correspond 
to minimum v. and the expression for 1/ ky is written and differentiated 
more simply. Thus, differentiating the right side of (11.4) with respect to 
0, and equating it to zero, we obtain 


b o - ® 
2 cos 294m cot > - sin 29 1m (2 + hy + h,cot? 4 = 0 


from which, after the left and right sides have been divided by cos 20m? 


we can determine Fim thus: 
d 
2 cot 3 
(11.8) tan 291m = taht! 
+ hy, + 400 5 
or 
h 
® 1 
; ) = — + 
(11.9) cot 2 1m tan 9 roe 


Let us now see if the obtained value of Gim corresponds to the 
minimum and not the maximum velocity. To do this, we use (11.5) to 
determine the directions to which an infinitely large value of the magnitude 
of the velocity corresponds. Evidently, these directions correspond to 
those values of 4 which make the denominator on the right side of 
(11.5) vanish. Equating this denominator to zero, we obtain 


6. Orein iO cot —heoa 6. (9 hen eee |) 
cos 0, sin 4, cot > cos 4\2 + h, + 1° 5 = 0, 
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! 


from which the desired values O46 and 04 ») can be represented as 
follows: 


cos 1. =-0 05. = ’ 
- - 
7 2+ hy + h, cot* > 
(11.10) - tan a = aa 3 
co 5 
hy 
- ta ai = ia 
1 Leo = ND + Gin & 


The physical meaning of the directions corresponding to the values of 
0, ., and 056 obtained is the following. Evidently, for 9, close to 6, 
the trajectories approach a vertical, rectilinear trajectory and for 05 
close to ¢]_ the trajectories approach a rectilinear trajectory passing 
through the points O, and Oo. That the angle 05. determines the 
direction of the line O,O9 can also be shown geomeirically by comparing 
the expression (11.10) with the relation which can be obtained from 


4 A009 (Fig. 8). 


wa) 


Comparing formulas (11.9) and (11.10), we can write 


t = us _ _ 1 a 
01 ~ 9 251m - Oo oo 281m 


or 


= 6 6 = 


ni 
=e lm 9 


(11.11) 9,, - 9% eee: 


1m 


Hence, it follows that the direction corresponding to the angle 0,,, is the 
direction of minimum initial velocity and it bisects the angles formed by 
the directions of the infinitely large initial velocities. 


The same result can also be expressed thus: The direction of 
minimum initial velocity which will guarantee incidence on the point S 
from the points O, and Oo bisects the angle formed at the point O, 
by the vertical and the direction at the point of incidence O,. 


The analysis of the family of identical-range elliptical trajectories is 
of great interest. 
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Using the property of an 
ellipse, which is expressed as the 
fact that its tangent bisects the 
angle between the radius vector 
drawn from one focus to the point 
of tangency and the continuation 
of the radius vector drawn from 
the other focus to the point of 
tangency, and taking account of 
the result obtained above, we can 
arrive at the conclusion that the 
second focus of an elliptical 
trajectory corresponding to 
minimum initial velocity lies on 
the line 0,05. The position of 
this focus, denoted by F in 
Fig. 8, is determined on the line 
0,09 by the following obvious 
Fig. 8 relation: 


> 
O 
— 
+ 
) 
© 
I 
> 
O 
ww) 
-}- 
Fr 
O 


or 


FO, - FO, = AQ, - AQ, = 


This relation yields a simple method of determining the position of the 
focus F onthe line 0,09 and of constructing elliptical trajectories 
corresponding to a minimum initial velocity by means of the points A, 
O; and Oo. 


Let us note yet another property of the elliptical trajectory under 
consideration, which is obtained easily by using the property, formulated 
above, of the tangent to the ellipse and by the fact that the focus F lies 
on the line 0,059. It can be formulated thus: The tangent, at the point 
Oo, to an elliptical trajectory corresponding to a minimum initial 
velocity bisects the angle formed at this point by the vertical and by the 
direction at the point O;- 


Using the last property, we can write (see Fig. 8) 


_  £0,0,0, 
ea ae 
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It follows from A4A0,0, that 


t a t rt ae 
LZ O,0,02 = + 6) + OF = &+ 26, 


As a result we obtain 


d T 
(11.12) 6,,, + [om + pa 


where 95, is the angle of incidence of the elliptical trajectory 
corresponding to the minimum initial velocity. 


Finally, substituting the obtained value for Gam into (11.4), we can 
obtain the value of the minimum initial velocity Vim: First, however, in 
order to simplify the calculations it is advantageous to use (11.9) to reduce 
(11.4) to the following form: 


1 ee 
(11.13a) ky = cot > (sin 20, + 2 cos*64 cot 264) 


after which we obtain 


} : 
Kim = cot (sin 294m + 200861, cot 2041) 


or 


2K am 


d 
(11.13b) Kim = tan ry tan G1m > Vim: = ry 
After the analysis of the properties of trajectories corresponding to 
minimum initial velocity, which we shall simply designate as minimum- 
velocity trajectories, let us return to the family of identical-range 
trajectories. 


It was established above that assignment of the value of the selected 
parameter O4 uniquely determined the trajectory corresponding to this 
value. However, it is perfectly clear that any other parameter, say, the 
magnitude of the initial velocity, can be chosen as the parameter of the 
family under consideration. In order to determine the direction 0, ofa 
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velocity which will guarantee passage of the trajectory through the 
prescribed points O, and Oo for a given value v, of its magnitude, it 
is necessary to solve relation (11.6) (obtained above between v, and 91) 
for 91. This is done most simply by returning to the above-mentioned 
dependence in the form (11.4). Dividing by cos, and multiplying by 
k, on the right and left sides of (11.4), we can write: 


ee 2k t = tan 6. + k 2+h, +h t” = 
or 
Qo | 2 : 7 f anboe Ve 


from which 


Z ® 2 422 ee 
tan 0, = k, cot > + k, cot E+ ky (2 + By + Boot EY 3 


or 


h 
® 2 ® 6 “I b 
(11.14) tan 0, = k, cot 9 cot 9 ky + 2k, tan (tan 2” oi ; tan 2° 


The written formulas show that for each value of kK, satisfying the 
condition 


h 
® ® 1 22 
rion $ (tan + 42) — tam? S > 


(11.15) kf + 2k 


there are two values of #, and, therefore, two trajectories corresponding 
to this value of ky and passing through the prescribed points O, and 

Og. One value of 4, and one trajectory passing through the prescribed 
points O, and Oy correspond to the value of Ky satistving the equation 


h 
2 2 2 J 22 
(11.16) ky + 2k, tan 9 (tan - + ain 5) - tan oo 0 
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Let us determine the roots of (11.16), which we shall denote by Kim: To 
do this let us use (11.9) and let us rewrite (11.16) thus: 


2 2 22. 
(11.17) Kim + 2k,_,tan 9 cot 204m ~- tan 5 0 


after which the expression for the desired roots has the form: 


rc) ro) ce) 
Kim = tan , cot 291m + /tan 9 cot 204m + tan 5° 


It is not difficult to see that the value corresponding to the minus 
sign is superfluous since the parameter ky cannot be negative. Hence, 
Kim has just the one value 


od 
= eee 2 
Kim = tan 9 (- cot 295m + / cot 294m + 1 . 


This latter expression is easily reduced to the form 


m 


® 
te: = tan tan O1m 


k 
which shows that the found value of Kim agrees with the value of Kim 
determined above which corresponds to the minimum-velocity trajectory. 
Hence, we arrive at the conclusion that there is a unique value of 6,, 
which equals 9 ,, as follows from (11.14), only for the value of kj,, 
corresponding to the minimum-velocity trajectory. Therefore, a single 
trajectory passing through the prescribed points O,; and Oo corresponds 
only to the value Kim and this single trajectory is the minimum-velocity 
trajectory. 


For each of all the other values of k, satisfying the condition 


(11.18) k, > ky 


to which the condition (11.15) written above reduces, as is easy to see, 
there are two values of 01; which we shall denote by 9,4 and 940 and 
shall define in conformity with (11.14) and (11.9) by the formulas: 
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Q = — ky + kK? + 2k eee 26 - age = 
tan 11 7 co 9 1 l 1 9 Im 9 
(11.19) 5 3 3 
= = 2 2 a 2 
tan O49 = cot a - Ky + 2k, tan 9 cot 295m tan a )° 


The dependence between the angles O44 and G40: expressed by 
(11.19), can be represented in simpler form. To do this, using (11.19), 
we write the expression for the tangent of the sum of the angles O44 and 
949. We obtain 


= 


2k, cot , 


il 


(11.20) tan(@,, + 95) 


® rc © 
Pas pai = =a 
1 - cot 9 (tan 3 2k,tan 9 cot 20, 


or 


from which there also results the desired simple dependence between the 
angles O44 and 950 in the form 


(11.21) 91, + 915 = 261m 
or 
(11.22) 94, - 9m = 8m - M49 


The latter relation shows graphically that two values of 6, 
correspond to each value of k, satisfying condition (11.18), namely: 
O44 and 940) which determine the directions of the initial velocities, 
symmetric with respect to the direction of the minimum initial velocity. 
Two trajectories passing through the prescribed points O; and Oo and 
corresponding to the very same value of the initial velocity correspond 
to these two symmetric directions. Hence, that trajectory which is 
located under the minimum-velocity trajectory is usually called the 
curved trajectory and that which is located above the minimum-velocity 
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trajectory is called the flat trajectory. Hence, the minimum-velocity 


trajectory is a sort of boundary trajectory separating the curved and flat 
trajectories. 


This latter property of the minimum-velocity trajectory permits the 
expectation that introduction of the parameter 


(11.23) Ae, = @ 


1” “1m 

in place of the parameter v, might simplify the original relation (11.4) 
for the family of identical-range trajectories. To verify this hypothesis, 
let us return to the relation (11.13), which is identical to (11.4), but more 
convenient for our purposes. Performing simple transformations with 
trigonometric functions, let us reduce it to the form: 


g Cos 2(9, ~ oie) + cos 26 


(11.24) + = cot = = 


after which it is possible to transform, without difficulty, from the 
parameter 0, tothe parameter 4@, and to write 


sin 20m S 


1 ~ cos 246, + cos 28 vanS 


(11.25) k 2 
1m 


or 


2 
2k COS ae 


(11.26) Ky = cos 248, + cos 294m 


Using the fact that the denominator of (11.26) can be represented by 
means of simple trigonometric transformations in the form: 


*In this very graphic definition for elliptical trajectories, we speak only of the 
sections of the trajectories included between the points O, and Oo and located 
on a sphere of radius ro. To apply this definition to the parabolic and 
hyperbolic trajectories, it is necessary to consider the location of sections 
containing the infinitely distant point under the minimum-velocity trajectory. 
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cos 2a + cos 28 = (costa - sin%p) + (cos2B - sin*q) 
= 2(cos*a - sin“g) = 2(cos2B - sin*a) 


let us rewrite (11.26) in the form: 


kjmcos76, 


(11.27) k, = 
. cos“ - sin2Aé@ 
lm 1 


Formula (11.27) shows that the expected simplification of the 
original dependence actually occurs. 


Sometimes the following relations will also be useful: 


| ar, mo 2 
Vv Nien kK, Kim sin Ae, 
Vim lm cos“6, - sin?A6, 
and 
k cos“6 1 - tan @ tan Aé 
1 1 7 1m 1 
ky m°08"94 mn 1 + tan Om tan AS, 


which are easily obtained from (11.27). 


In conclusion, let us examine how the doubled sectorial velocity o 
changes during the transition from one trajectory of the tamily to «nother 
To do this, we put (4.9) for © in the following form: 


e a 2 
(11.29) o = rjV,cos a, = 27 or 4k, cos o, 


It follows from (11.29) that the change in the doubled sectorial velocity 
from trajectory to trajectory is determined completely by the quantity 
kycos0 1, since the remaining quantities remain constant for all the 

trajectories of the family. On the basis of (11.27) we can write 


ue 2 
Kym) C08%81,,,C08 0, 


lm 
(11.30) k,cos*0, = 

1 1 9 ; 

cos“@é - sin“dé 


lm 1 
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and, having made simple trigonometric trunsformations, we can represent 
the expression for k,cos“0, in the following form: 


ky 


1.8) (es) 
(1 - tan osm) + 2 tan Gam tan a, 
In (11.31), the expression 1 - tan*41,, is always greater than zero; 


9) : 3 ‘ 
hence, k cosUy diminishes as 0, increases and, therefore, the doubled 
sectorial velocity o@ diminishes. 


The expression 


“le 


can be written for the flight time T, where S is the area swept over by 
the radius vector during the motion of the point under consideration over 
the trajectory of the family; o is the doubled sectorial velocity of the 
same point, which is constant for each trajectory of the family. 


It is clear from geometric considerations that the area S becomes 
large as aT increases. On the other hand, we explained somewhere 
above that the doubled sectorial velocity 9 decreases as 9, increases. 
Therefore, the conclusion can be made that the flight time over the 
trajectories of the family increases monotonically as 94 grows. 


We thus terminate the analysis of the properties of identical-range 
trajectories. Let us just examine the following important property of a 
minimum-velocity trajectory, which we have not touched upon. The 
minimum-velocity trajectory of a family of identical-range trajectories is 
at the same time a maximum-range trajectory for a family of identical- 
velocity trajectories. Starting from this property, we introduce completely 
identical notations for the parameters of both trajectories since the 
relations obtained for the parameters of one trajectory are also valid for 
the parameters of the other because of this property. 


The identity between the above-mentioned trajectories is proved 
most simply if, when analyzing the family of identical-initial-velocity 
trajectories, the magnitude of the initial velocity is taken equal to the 
magnitude of the minimum initial velocity of the family of identical-range 
trajectories. Evidently one trajectory among the family of identical- 
initial-velocity trajectories will absolutely pass through the point Oo 
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prescribed for the family of identical-range trajectories. This will be 
the trajectory which coincides with the minimum-velocity trajectory. Let 
us prove that this trajectory will be the maximum-range trajectory for 
the family under consideration. 


Actually, if this were not so, then among the family of trajectories 
there would have to be some whose range exceeded the range between the 
points O, and Ops. In this case there would certainly have to be one 
such trajectory which would pass through the point Oo; but this contradicts 
the original assumption on the minimality of the initial velocity taken for 
the prescribed points O, and Oo, since only one trajectory passing 
through these points exists for the minimum initial velocity. Therefore, 
the minimum-velocity trajectory of the family of identical-range trajec- 
tories is at the same time the maximum-range trajectory for the family 
of identical-initial-velocity trajectories. Hence, all the relations 
obtained for the parameters of one trajectory are valid for the parameters 
of the other as well. 


12. FUNDAMENTAL RELATIONS BETWEEN 
PARAMETERS OF THE MAXIMUM-RANGE AND 
MINIMUM -VELOCITY TRAJECTORIES 


Because of the great practical interest represented by the maximum- 
range trajectory of the family of identical-initial-velocity trajectories 
and the minimum-velocity trajectory of the family of identical-range 
trajectories, it is expedient to present in a separate section, the relations 
between the parameters of these trajectories, obtained in sections 10 and 
11. This section is in the form of a summary of the formulas intended 
for the solution of various practical problems. 


If the parameters 
(12.1) ry, To, ® h, See 


of the mentioned trajectories are known. then it is convenient to define 
their remaining parameters by means of the following formulas: 
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rt 
2 cot ry 
(12.2) tan 284m = 


™ rea 
2+ h, + h,cot 5 


My *9 Yim’? Kim oT 


are known, then to determine the remaining parameters it is convenient 
to use 


(12.0) tan ae = 


(12.6) tan— = k, cot @ 
m 


13. DEVIATIONS IN THE ANGULAR FLIGHT RANGE 
CAUSED BY INITIAL DEVIATIONS 


When analyzing motions along Keplerian trajectories it is often useful 
to know the dependence-relating deviutions in the initial-motion conditions, 
which we shall simply call initial deviations for brevity, to the deviations 
in the angular flight range which they cause. 


This dependence is given implicitly by (7.18), which represents ® 
as a function of the quantities Vy 94, ry: 


(13.1) & = P(v,, 9, F)) 
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Explicitly, it is usually represented by the Taylor formula* 


nN 
ae k 
(13.2) 5 = one + Ray 


where 


k 
i ep a. age . 


Rn+1 is the remainder term of the formula; 5v,, 66 D dr, are initial 
deviations in v, @ and r; 6®@ is the deviation in the angular range due 
to the initial deviations OV1) 644, Or, - 


Considering the initial deviations 6v,, 601, dr, to be sufficiently 
small, one usually is limited to an approximate expression for the desired 
dependence; terms higher than first and second order infinitesimals are 
neglected in the Taylor formula. 


1. Expression for 616. The expression for 61 can be written in 
the following form:** 


(13.4) 646 = @ dv, + $ 66, + %, Sr, 
1 1 1 
where 
_ v® =_ JY® _ vd 
(35) & =“ y,,% 2-2, 3 -2: 
: Ov, _ _ ory 
(13.6) ov, = vy 66, = 64, ; ory as 


where the angles are: measured in radians. 


*Here and everywhere henceforth the symbolic form ‘of writing the Taylor 
formula, which is convenient for its simplicity and clarity, is used. 


**As before, we shall denote a non-dimensional quantity by a bar above the 
letter. The bar will also be set over the angular quantities for symmetry in 
the notation. 
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To determine the quantities in (13.5) it is necessary to find the 
partial derivatives of © with respect to the parameters Vy: 0, ry: 
One of these derivatives was determined in section 10, namely, the 
derivative J® U05, by differentiation of (7.16). Let us find the remaining 
derivatives in the same manner. Let q, conditionally, denote the 
variable with respect to which we shall carry out the differentiation. 
Then we can write 


® = 
\ ple? 
— l tan 5 tan O, ab - a o OX, - ’ Oo tan a, oh, 
; 9 ® oq = 2 aq ot 2 oq aq - 
cos ry 


In order to obtain the desired derivatives by using the above formula, 
let us write the expression for the partial derivatives on the right side. 
We shall here assume q to be not only the parameters vj, O41: eT but 
also the parameters 7) and Lo, since the deviations in range caused by 
deviations in the parameters T™) and ro are also of considerable 
interest and can be obtained in passing together with the desired devia- 
tions. Thus, putting 


(13.8) q = Vz G4, ry, ™: Yo 


we obtain: 
= x 23 
pits at eek ae a: Oe es | 
: Ova rele OT y , 1 Yo OL re 
oO tan o, 0 tan o, 0 tan a 0 tan a, 0 tan a, i 
(13.10) —— = ooo Pm = oF 
0 
OX, 9 1 OMY 1 es ee: 1 
7 or. tT. ° 004 — 
oa k,cos*0, “1 “0 k,cos“#, 0 1 k,cos0, 
(13.11) 
eh, ee te 
oro ory 2 r k,cos a, 
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Substituting (13.9), (13.10) and (13.11) into (13.7) and using (8.8), let 
us write the desired expressions as follows: 


peiee 2 
ab sin 2 14 
na = 
1 2 
k cos 0, 1 
ab sin & @\ ( ad 
930. = Y1 a (&, - tan aT tan | = v1 (leo a tan 0) 2 
1 k,cos“@ 1 
1 1 
sine 
96 1 2 ee es | 
(isis e = Vila a oe ee 
1 2 k.r4cos“6 i 1 ro 
11 1 
re eed 
OOS ui ee ‘1 06 26 x 
— —_— 1 —_— — ———_——_—- ; ——— ~~ me 
2 
1 2 Ei. 7 5 
Q, = h, = 
tan 6, + h es : -) 
1 1225 


As a result, the expression for the deviation of the angular range ble 
caused by deviations in the parameters (13.8) can be written as follows: 


2G Oi a ér, + @, 67, + $, oro 


(13.13) 646 = 6 dy 
= 1 o 9 2 


1 


where 
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_ 3 2 sin” — 
Vy avy “1 ~ F1 
: 1 k,cos 0, 
} o° (. t tan 6 S 
6 = Wa = co ~ tan 
1 00, 1 1 vy | 
of = | 
(13.14) <@ =— r= @ + ~. — 
vr) ory 1 9 vy 1 Po 
3 - & __.1ig 3 22, j lL 
r. - 7 roa 
tand:, whee ry Lr, 
1 1o° 9 
7 Ov _ : or, 7 Oro = OT 
(13.15) ov, = ——; 66, = O04 (Ore = Or, SS 5 01g. = 


2. Expression for 56, The expression for 6° can be written as 
follows: 


Oe a “9. = ~2. -2 . =x - 
(13.16) 6° =8  6v2+%, 60%°+ 5  Sr7+H | bv, 68 
® 6 + ® 5646 
V4"y elec | 0,14 1°Ty 
where 
7 L 2 Ae 5 La’e eal a2 a° 
ee ’ oR : = 1 
Vi 2 1 oye 0,9) 2 962 rir, 2 or’ 
1 1 1 
(13.17) 
5 9” - 9° 9° @ 
v19, 1 0v,06, © vyry Vl Ovyery © 9404 00, or) 


and dvy, 60,, 6r, are defined, as before, by (13.6). 
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To determine the quantities in (13.17) it is necessary to find the 
second partial derivatives of ® with respect to the parameters rj, Vy: 
6,. This is done most simply by differentiating the first derivatives 
(13.12) obtained above. Hence, as before, we determine the second 
partial derivatives with respect to all the parameters of (13.8) in passing. 


Returning to (13.12), we note that the majority of the first derivatives 
are expressed in terms of the derivative d¢/ OV1- 


Hence, by determining all the partial derivatives of this derivative, 
we facilitate the determination of the remaining second derivatives. By 
denoting by q_ the variable with respect to which the differentiation is 
carried out conditionally, we can write, starting from (13.12) 


(13.18) ort. e a “1 oes 2, jeg ate 
0 14 Ov Py aa sin? = AY 11 Aq kyvy 

+ cos“6, ba ; 

cos“ 6 


Substituting the expression for ay, /8q here, which is easily 
represented in the form 


6) 
P P hy 9 é O 1 0 tan 6} 
(13.19) oo 1 Se ay he ee 
2 sin” — a oq oq 
we obtain 
2 oh 
o ® 
(13.20) > = + + cot 2 ae 
14 2 sin* > 1 oq 
89 tan @ b dh, 
+ sn 26, ?4) aq ~ 4&4 cot 5 <i 


1 AK Vy) “Vp 
Kivy OY 


el 


Using formula (13.20), we can write the following derivatives without 
special difficulties: 


2 gh D 
Oo“ @ & 
(13.21) = cot > + a (22) = 3 Od 


avi asin2=/ \V1/ My OVy 
2 
9° 6 & "5, \5e 90 % \oe 
(13.22) aoe, cot — + av. 86 + /2tan 6, - —- 
vl 2 sin? — 1 1 cos"6, ovy 
; 7 
is x ae & 75, \56 90 %  o6\oe 
(13.23) 5—— = {cot > + ae ~{ = + — cot — )— 
1 osine—/% 8%, \T1 2 8 2/ dvy 
2 oh 
1 0 2 sin” > 1 ™ ™) ovy 
2 yh . 
ee: & 111 \ae ae 1 & 9@ 
(13.25) = (cor $ + “Ha, a S49 ct— = 
OV Or5 2 2 sin?= Ov, or, Lo 1 2 OV, 
_ <i ab ae 
2 r 


Now it is convenient to return to the derivative 8¢/ or, and to 
differentiate it with respect to all the parameters of (13.8) except Vy: 
Here, as in the preceding case, it is convenient at the beginning to write 
the general expression for the derivative of 0¢/ or, with respect to q. 
Using (13.19) this expression can be represented as 


2 


~ 
a“ ® os | Naas | 1 2 
=e - —= 
(13.26) ar. 8q ,2 oq 2 84 3 oq 
y) 2 2 
h h 6 
Es ee fae le | Ve One 
To 1 aie 2 oq oq ry oq ry 
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Using (13.26), the following derivatives can be written without 
special difficulty: 


ae > 1 96 26 4 ae 
or ral) 9 Sine 2 OF 5 Q 1 és 26 Yo 
D 1 
on “5, \ oe ae 
® ] dr, 96 
2 sin” — aca | 
9 
?. + cot t 2% \ 
cos~6 2 00) OF 5 
1 
é “01 a6 a6 1 % _@\a@ 
22) ee - Scot >} a 
%5, \5e 96 
= |cot— + B aa or 
2sine?—/ 72 1 
2 
(2 ee Loy & FE 
Ty lo 2 2ry 2 og, Lo 
26 Pb, 9g 9 
(he:22), ae. = b or. dr, 
2°"0 2 sin*—> Rg: 206 
gh 
11 96 ao 9b 96 
. core or. On 20 or 
2 sin” 2 0 To 0 
(13.30) ve = as i ) oe cot - 2) 28 
or? 2 sin2® \p 2 2° ro} ar, 
2 
ai Pp ‘ hy ad ) @ Vb Od 
Cc 7s ) "Yas aa + basa CO as ie ae 
os sin sa (ar ° i) OP» 


Furthermore, let us differentiate the derivative ab/ om with respect 
to all the parameters of (13.8) except Vv, and r,- As a result we obtain 


2 ie 9 
(13.31) ~ = re te ee 
yh P 
b 
= (o§ + Ts os Be + (een, 
(13.32) = = = ee 
19 ri 279 OV4 Ory 
& gh , v) 
= cot a ape ol OP 
c 2 sin” = Oy Oy \Ty To 2] OM 
RE ge DO Oo ies (22-) 
ane 27 8V1 9% ane Ovy 2 6 sin” 2 aT 


After this, let us write the expression for the derivative 9®/ dr, in the 
form: 


Ge, pees: Sy eis e, 
(13.34) = 


and let us differentiate it with respect to a, and r,;. Asa result we 
will have 


es ao  ‘y of6 "2 ate 
ér 00, 2ryj dv, 08, ry Or 98, 
fot #5, \ae ae * 90 
7 ar, 00, — 8 
2 Since a eae Cc s°0, "1 
2 
Vy rT ns a el 
+ — tan 0, — --— cot- = 


@ Vv 26 Yo 9%6 1 oe 
(13.35b) 22 = pl ee Ge ee et eee 
2r, oOv.,or r Or,or, ry or, 


2 
or 1° 1 1 
1 
fe. Yb, \o@ 90 (2 oe #08 20, 
= {cot — —— rns — — |— 
Seine = or, Or, ry Yo 2 ory r4To 


We determine the last remaining derivative by differentiating 
O6/ 86, with respect to 0, . After differentiation we obtain 


9° o Viki ae ae 


2 
a3.se) 22% =, (. re tan ¢,) 
ris fe ta | = 1 6. & 90 
06 2 Nea pine 
2 
‘1 ae 
cos*6, ovs 


Substituting (13.22) here, we rewrite the expression for the desired 


derivative as: 


2 gh 

Qo ® i) 11 \0¢ co) o® 

(13.37) cot 5 + 9 & | a8 tv, (cot i tan | * \ 
si S 1 1 


2 
06, 
My 
1 ® o® 
{ra (seo! a tan | se | 


+ {(2tan 6, - 
cos“é 
k 
_f_“r'1 ae \ae “1 ab 
. a 


9 


o 


in the parenthesis by its expression 


Replacing the derivative 0¢/8v 
O® 06) and making 


and the expression in the braces by bh 
other obvious transformations, Ict us write the expression for the last 


derivative in the following form: 


e derivative 


98 


2¢. Vv 
+ 2tan 0, - ——— be) dP 


06 
cos", 1 cos 1 1 


Hence, an expression can be written for the deviation 6*® in the 
angular range caused by deviations in all the parameters in (13.8), as 
follows: 


2 -2 = 2. -2 <= - > 
(13.39) 6¢6 = 6 dye + o O65 + @ orp + ® , 5v456, 
viv, 1 11 ae Fa | v494 


+ 
oO 
mS 
< 
pany 
& 
a) 
—_ 
+ 
Hl 
(or) 
® 
—_ 
en) 
ar) 
_ 
af 
1 
eu 
pay 


159) Tory ory 
where 
ve 12% 2.1 gs 3 -33 
( ; ) V1%1 2 avi 1 2 2 Va V1 2 V1 


2 P 
= 1 oo en ee 1 os 
gy Pe 6. -———- | 
(13.41) ° 68 => 7 5 2g , + (an 1 5 9 
1 
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2 
= 120° 1 
== ri > =5g8 6 + g(lth 
(13.42) 2. 2% 2 2%2°r % ¢ 
161 ory 1 
® 
= 1 + g(1 + h,) cot =) oo 
2 
cas 9° 
(13.43) © =v = 6’ +/2tan6, - 
V4 94 1 v1 99, 2 vy 0, 1 
2 
- 9°@ Re : ie 
(13.44) @ aViten = ee [1+ qa + Boot $a 
V4T] 11 Ov, or, 2 V4 ry 1 2 
= 26 Se * 5 
(13.45) © =rY 2006 
0404 1 06, or, a4 Ty cos?6, ry 
= db — aa 
z 9, (1 + h,)cot 5 ®, + tan a, e 
1 1 
2 12086 1 -. = 
(13.46) @ = a eer a fp? oo 
“00 am 0 “oO 
2 
5 ® = b — 
(13.47) @® = ae — = + 8 & - (1 + cot > ® 
Fore ars Yo To Lo 
2 
= 9“ = © _ 
® = = _ ® 
(13.48) TV TP an aT 5 a to cot — ® 
es 
2sin?-— "0 2 
(13.49) & To a*@ oe ae 
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= od = 
(13.50) © = lov. == GP OCG 1+h cots & 
gh 
sa = 5 5 
2sin*-— "2 4 
2 
2 
= 9“ a ? 
(13.51) 3 = ~—- = 9 @, +(2tano, -—-)\ 3 
1,0 0 d7.96 2° 1, @ | 
071 Maal 0 74 cos“6, T 
2 Y 
| * a“ 6 = = 1 6 9b \ - 
(13.52) i = Lee ee ee rc) rc) = + cot ——\ @ 
r.9 2 or.00d 2 g 
2"1 20% "2 "1 \ cosa, 2 00;} Yo 
(13.53) 8 zor o"® yd 8 fa + (1 + hy) tls 
; —— el z g, (1 + hy)cot — 
ToXy 01 amy Ory 2 T) Vy 1 1 2} % 
2 
_ oO“ 
(13.54) & = ror} 
ro0y 2 or,or, 
ee. gis ce mo | ®@ — _ 
= 6 © - St a 
Yo ry r + 2 29, (1 + h,)cot 9 9 cot 9 f| - 
oh 
d 
(13.55) % = cot> + ee a ;  g = ——+—_——_ : 
2 sin ry tan O4 + hycot > 


and Ovi» 68,, or, Om, or, are defined by (13.15). 


In certain cases the deviations in angular range caused by the 
initial deflections are conveniently represented in the form 


2 = = 1 y) = 1 es = 1 aoe 
(13.56) 6° P@ © Go + *00,° $08, + “or, bor, 


to which (13.16) is easily reduced if the substitution 


ae 1 1 = ae 23 = 
(13.57) 6v, = (° b- &, 60, - sr, | 


is made on the basis of (13.4). The following expressions are here 
obtained for the coefficients in front of the products of the deviations: 


= 1 3 1 
(13.58) 5,75 % ~ 9 5 
1 


5 

1 

= 2tan 0, - 9 (1 + tan%@,) + 3 = 
V 


(13.59) * 36, = 
1 
: 2 9 2 
= 3k, cot 2 7 tan eo P, (1 + tan 95) 
d 
= = ® aT 
(13.60) "er, = -|1 + (1 + hy)eot | + 3 ca 
1 
ere ay ce 2 
= 3 + (1 + D(a - cot | 
V 
1 
*, 
(13.61) & = - aye + tan“@ )® - tan 6, © = 1 
6,9, y) l’"v 1 $86 2 ® 
1 1 V 
1 
1 3 ®\ . 1 
= (—t 0. -— — =. 2 
(2 fan 6, - 5 k,cot a pil + tan“6,) ey 
%, 
(13.62) Parry = sin (1 + tan“6,) - 2tan 6, ®,. - 3 = ?. 
V 


1 


2 : @\ _ 
Y (1 + tan 9,)sin @ + (san a, ~ 3k, cot 3) oF 
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bo [co 


(13.63) on = @(1 + h,) 


14. DEVIATIONS IN FLIGHT TIME CAUSED BY 
INITIAL DEVIATIONS 


In analyzing motions along Keplerian trajectories, we often use a 
dependence relating the deviations in the initial motion conditions to the 
flight-time deviations they cause. Applying the Taylor formula, as in the 
preceding case, to the function 


(14.1) T = T(vy, 94; ry) 
we can obtain this dependence in the following form: 


k 
(14.2) 6T = 2 ST+R, 


where 


1 9 9 8\k 
(14.3) o8T = ia ("1 By, + 66, 3e * or, =) T 


Rn+1 is the remainder term; 6ov,, 68,, 6x, the initial deviations in the 
quantities v, r, 9; 6T the flight-time deviation caused by the deviations 
OV1> 66,, or, - 


Considering the initial deviations to be sufficiently small, it is 
possible to find an approximate expression of the desired dependence by 
neglecting terms higher than the first or second orders in Taylor's 
formula. 
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The expression for 617 can be written as follows: 


a = = 7 ee - 
(14.4) 6°T = Ty ov, + Ty 06, + a or, 
1 1 1 

where 

== OT a OT oT 
(14.5) T A Vie Dip: ee =. oe 

vy avy 1 1 00, ry or, 1 

ov, : or, 

(14.6) Ov, = Vv, ; 06, = 68, ; ory = 7, ; 


In order to determine the quantities in (14.5) it is necessary to find 
the partial derivatives of T with respect to the parameters Ty> o> V1: 
Let us determine these partial derivatives for the case of point motion 
along elliptical trajectories when the function (14.1) can be represented 
in the form of the relations 


T TT 
- = «1, |(; = Ry * ”1) “(5 - By + "2)| 


08S pee 
p T™ 2(1 - ky) 
1 - 2k, é 
B, = arcsin (- — = By =2| 
(14.7) 
1 - 2k, 
Bo = arcsin (- 5 = Bo =1) 
. vir, k, + h, 
¥y = € COS Bp, ; k, = : ky = 
on 1 +h, 
ie. 4! 
A 1 2 
Yg = € COS Bo ; hy ease aria e% = (1 - 2k,)"cos” 0, + sin*6, 
2 


presented in section 9 in formulas (9.59) - (9.61). 
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For our purposes it is convenient to represent the first of the above 
relations in the form: 


(47a) T= CH™T [ry + (-1)7,] 


where 


> 


: TT 
(14.8) No tT V1 7 Py: T = 5 + Y2 - Bo 


Values of the numbers m and n must be selected as follows: 
1) If the points O, and Oo between which the flight time is 
measured lie on different sides of the peak of the trajectory O.,, then it 
F 
should be assumed that 
(14.9) m=0; n=Q 


2) If the points O, and Op, lie on the same side of the peak of the 
trajectory and ro < ry, then it should be assumed that 


(14.10) m=1; n=1 


3) If the points O, and QO, lie on the same side of the peak of the 
trajectory and ro > rj, then it should be assumed that 


(14.11) m=0; n=1 


We find the derivatives of interest to us by differentiating (14.7). 
Hence, as in the preceding cases, let us determine the partial derivatives 
with respect to all the parameters. Denoting the variable with respect to 
which the differentiation is to be carried out by q, we can write 


OT OT OT 
OT m 1 ge 2 


It is easy to note that 


tiny Pe [2 2 
(14.13) SO = “pl 2a aq 2m, 8q 
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Before determining the derivative 9a/éq it is convenient to write 
the derivatives of k, and ky: 


ok OV or On 
1 1 1 1 0 
(14.14) ae = (2 ee ee ee eee ca ee =) 


oq Vy od ry ™; 94 
- 10 
oq 1 + h, v4 oq | oq "9 oq a) oq 


Differentiating the expression for a, we obtain 


or ok 

da 1 1 1 1 

; — ee a ial eee iaatiemiatetimeaeaial + | iiiniintiieail et ee 
(14.16) (2 ; i - a 


Substituting (14.14) into (14.16), we can rewrite it 


aq rl - k,) aq ~ v,(@ - k,) 84 — mL - k,) aq 


Using (14.17), we can represent (14.13) as follows: 


p 
(14.18) —— = u : 


oT or 3k ov. 1+ 2k, on 
yy | Seaee n OO ORS “| 
oq Plar,(1 - kj) 94) vy(L- k,) oq 2m(1 - ky) oq 


from which 


OT 
ae 3 a T i = 0: Mar = es A T 
1 v, 1-k, P 06, aT 4 2m,(1 - k,) P 
(14.19) oT, ; : . aT . 
or, 2r, (1 - k,) p OL 


Differentiating the expressions for Ty) and Y¥z> we obtain: 
OT ay ap 


1 1 
14.20) —=— 
( oq oq oq 
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oy op 
1 _ oe 1 
(14.21) a Gd cos B, - esinp, Bq” 


after which, substituting (14.21) into (14.20), we can write: 


9 
it ta ah «A * ety 
dq aq 1 


1 
1) aq 


Differentiating the equality 


1 - 2k 


. = 1 
(14.23) sin By = z 
we can write 
B, SMP, 9¢ go 


eee) ‘€q. se cos By oq 7 ecosfB, 


Substituting (14.24) into (14.22), we arrive at the expression 


a 
"% 


1 de on 
(14.25) “a = @ cos By (sin 8; + e) oa + 2(1 + e sin B;) q. ; 


or taking into account that 


(14.26) Y= € cos By ; 1 + e sin By = 2(1 - k,) ; 


at the following expression: 


OT Ok 
i eee a ae be 7 —_t 
(14.27) oq Y4 c B, + e) aq + 4(1 k,) = . 


We obtain a similar expression for the derivative of To analogously: 


OT 1 de ok, 
(14.28) eq = Lom By + e) 8q + 4(1 - k,) “8q_ 
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Differentiating the expression 


(14.29) e* 


let us write 


de 4k, (1 - Kk) 08 ; 
OE ch. eg ; 9. —_——_- — - G 
(14.30) bq sin 0, cos 4% = aq 2 cos 4, 
Let us use the notation 
‘ 1 - 2k, 
(14.31) 7 =7, + (-1) T. 3 b, = sin By = 2 
1 - 2k 
b, = sin Bo = ee 
2 2 e 
then 
OT OT b +e b 
OT 1 n 2 oe| 1 n 
14.32) —=— + (-1) —=-> + (-1 
Sai ale ee a 
1 - k. 0k - 
+4 in aes 2 Ae 
or 
ok n ok 
a4a3) aaa - ky] t —. 1 2 
de D, +e " b, +e 
an + (-1) 
oa} 3 Yo 
since 
1 - k 
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(l - 2k, )"cos24, + sin0, 


1 - 2k 


e 


Using (14.33), (14.30), (14.14) and (14.15), we can write 


oT _ 2 OT 1 


k n 
OT 1 = 
or 4 7 Bs y= 40 - ko + “3 


or r ae 
Y - 2 
1 1 1 (+h) %2 
(14.35) 
b. +e 
OT ») 1 
sd OU (“Kn + DS); Da = 2 cos” 6 of 
OT» 0 1°3 3 3 1°11 1 
b +e 
+ cap 22 
Yo 
nt1 (1 - k.)* 
OT, 2h (ER fe ee, 
orp *2 (1+ 8)? Y2 


Then substituting (14.35) and (14.19) into (14.12), we obtain the desired 
derivatives in the following form: 


k 
OT 1 m 1 
ov, E t= oe ee *) % 
1 1 1 
OT m 1 
oT 3 m 1 
(14.36) or, E- - Ke) T + (-1) Ty(%q - 2) ry 


oT 1 + 2ky m 1 
<—— = . say T+ (-1) Thl-Ky¥3 + *) 4 


2 
oT yytnt coma Ue 
oro (1 + hy)? Yo Tp 
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where 


k, = 2m, 
; : k +h 
1 e ° 2 ie 2: 
1 
r - r 
- = 12 e = aa = : 
hy = To 5 ea/l 5a he oF “1 
r 
; a 
a=o, -n): Ip =4/7_ 
21 - k,)’ ~P "0 
| | - . T Ducsee 
B, = aresinb, ; By = aresin by Ca ee 
(14.37) 
m " - 
T= U™T|r + M7]: m= E+ ry - By 
To =o + Yq - Po 
sec. 
Ve SEL er ky nly. 
1 @ +h)” °2 
k 1)" 1 
Yq = 4(1 | en 8 
1 @+hy)* "2 
fo +e n Py + © 
b, = 2b, k cos a, yy td) Vn 
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As already mentioned, in certain cases the deviations in flight time 
caused by initial deviations are conveniently represented thus: 


(14.38) o'r = 7 646 + fF, 60, + BP. or 
9 6 1 r 1 
1 1 

to which (14.4) is easily reduced if the substitution (13.57) is made 
therein. The following expressions are hence obtained for the coefficients 


Ts Digs Be, 
a a 
T ® ® 
= “fe mit 3 7 . 1 
(14.39) Ts a 4 ae = T - TY =: f; at. T, = 
Vv 1 1 1 v 1 1 1 v 
1 1 1 
where 


Tr ere 
Vi 1- k 


m 


= m 1 
(14.40) Ty (-1) T,tan 9, ( + He 


a 3 m 

= ——$—=+« + (-1) T - b 
and the rest of the quantities in (14.39) and (14.40) are defined by formulas 
(13.14) and (14.37). 
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10. 
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